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Abstract 


I Closed-form  solutions  based  on  acoustic  theory  have  recently  been 

obtained  to  describe  the  wave  motion  in  both  low-pressure-ratio  rectangular 
(constant  area)  and  pyramidal  shock  tubes  which  utilize  different  driver 
and  channel  gases.  These  new  solutions  are  in  excellent  agreement  with  ex- 
perimental data.  This  work  should  be  of  interest  to  researchers  who  are 
using  shock  tubes  or  similar  devices  to  produce  impulse  noise,  in  particular 
the  simulated  sonic  boom,  in  order  to  facilitate  studies  of  the  effects  of 
inqpulse  sound  on  humans,  animals  and  structures.  Furthermore,  this  work 
is  relevant  to  the  understanding  of  the  wave  motion  produced  by  weak  planar 
and  spherical  explosions  of  finite  size. 

I 


Table  of  Contents 


Acknowledgement  s 
Abstract 

Table  of  Contents 
List  of  Symbols 

1.  INTRODUCTION 

2.  RECTANGULAR  SHOCK  TUBE 

3.  PYRAMIDAL  SHOCK  TUBE 

3*1  General  Solution 
3.2  Special  Solution 

4.  DISCUSSIONS  AND  CONCLUSIONS 

5.  REFERENCES 
TABLES 
FIGURES 
APPENDICES 


Page 

ii 

iii 

iv 

v 

1 

5 

15 

15 

30 

33 

37 


iv 


List  of  Symbols 


®1 

a2 

A 

^(t)) 

gid) 

\(3) 

Hit} 

m 

o 

M 

Ml 

P 

Pi 

p2 

Ap 

Api 

Ap2 

Apo 

r 

r 

o 


speed  of  sound 

sound  speed  of  the  channel  gas 
sound  speed  of  the  driver  gas 

cross-sectional  area  of  a rectangular  shock  tube 

1 n velocity  potential  function  of  the  driver  gas 
(inward  moving  waves) 

til 

i velocity  potential  function  of  the  channel  gas 
(outward  moving  waves)' 

til 

1 n velocity  potential  function  of  the  driver  gas 
(outward  moving  waves) 

Heaviside  unit  step  function;  H(t)  = 0 if  t < 0, 

H{t)  = 1 if  t > 0. 

mass  of  the  driver  gas 

molecular  weight 

molecular  weight  of  the  channel  gas 
molecular  weight  of  the  driver  gas 
absolute  pressure 

absolute  pressure  of  the  channel  gas 

absolute  pressure  of  the  driver  gas 

overpressure,  excess  pressure  or  pressure  difference 

overpressure  of  the  channel  gas 

overpressure  of  the  driver  gas 

pressure  difference  across  the  diaphragm  in  a shock  tube 
radial  distance 

radial  location  of  the  diaphragm,  length  of  the  pyramidal 
driver 

initial  radial  location  of  a fluid  particle  in  the  channel 
of  a pyramidal  shock  tube 


Vf 


v 


At 


fluid-particle  displacement  (r-r#) 
contact- sur face  displacement  (r-rQ) 
gas  constant 

universal  gas  constant  (Rq  = MR) 
time 

temperature 

temperature  of  the  channel  gas 
temperature  of  the  driver  gas 
fluid  velocity 
velocity  of  the  channel  :$as 
velocity  of  the  driver  gas 
volume  of  the  driver  gas 
initial  volume  of  the  driver  gas 

distance  measured  from  the  closed  end  of  the  rectangular 
shock  tube 

location  of  the  diaphragm  in  a rectangular  shock  tube, 
driver  length 

initial  location  of  a fluid  particle  in  the  channel 
of  a rectangular  shock  tube 


e (*,  t - x)/xo  or  (ag  t - r)/ro 

7 ratio  of  the  specific  heats 

7X  specific  heats  ratio  of  the  channel  gas 

7£  specific  heats  ratio  of  the  driver  gas 

i)  (ag  t + x)/xQ  or  (ag  t + r)/rQ 

*■'  (Vl  + a2P2)/ai(P2  " Pl) 


A wave  length 


I (axt  - x)/xq  or  (a^t  - r + 

p density 

p density  of  the  channel  gas 

P2  density  of  the  driver  gas 

<D  velocity  potential 

4>1  velocity  potential  of  the  channel  gas 

<t>  velocity  potential  of  the  driver  gas 

2 

factorial  function  (0 1 = l) 

(n)  binomial  coefficients  [m'./n! (m  - n) ] 


1.  INTRODUCTION 

Noise  pollution  is  becoming  more  significant  in  recent  times  and 
is  contributing  to  the  stresses  affecting  man.  Ongoing  studies  into  the 
effects  of  noise  on  humans  (e.g.,  see  Kyter's  book  - Ref.  l) , animals  and 
even  structures  axe  proceeding,  fortunately,  at  a corresponding  accelera- 
ted pace.  One  important  and  active  area  of  current  noise  research  i^  the 
study  of  the  effects  of  impulse  noise.  A particularly  important  impulse 
sound  is  the  sonic  boom,  whose  impact  on  society  is  being  assessed  before 
supersonic  transport  (SST)  aircraft  such  as  the  Anglo-French  Concorde 
and  Soviet  TU-144  are  introduced  into  extensive  commercial  service. 

At  the  University  of  Toronto,  Institute  for  Aerospace  Studies 
(UTIAS),  three  different  but  conplementary  sonic-boom  simulators  have 
been  developed  to  help  assess  current  societal  problems  associated  with 
the  sonic  boom.  The  first  of  two  major  facilities,  a loudspeaker-driver 
booth  (Refs.  2 and  3) , can  easily  accommodate  one  human  subject  or  a few 
small  caged  animals  in  its  solidly  built  and  sealed  chamber  (volume  of 
2ml)  to  facilitate  studies  of  human  and  animal  response  to  a simulated 
full-scale  sonic  boom.  The  second  major  facility  (Refs.  2,  3,  4,  5, 
and  6) , a travelling-wave  horn  in  the  form  of  a horizontal  concrete 
pyramid  (25  m long,  3-m-square  base),  has  at  its  apex  a specially  designed 
valve  which  regulates  the  air  discharge  from  a large  reservoir  into  the 
pyramidal  horn.  This  controlled  discharge  of  air  generates  in  the  horn 
interior  a simulated  full-scale  sonic  boom  for  human,  animal  or  struc- 
tural response  investigations.  Alternatively  a shock-tube  driver  can 
be  installed  at  the  horn  apex  to  produce  a simulated  short- duration  sonic 
boom  to  facilitate  certain  response  tests.  The  third  sonic-boom  simu- 
lator (Refs.  3 and  7),  a portable  shock  tube  (ll  kgm. , 1 m long)  having 
a constant-area  driver  and  an  exponential  horn,  can  be  easily  transpor- 
ted and  operated  by  one  person  to  conduct  wildlife  field  tests.  A 
simulated  short- duration  sonic  boom  can  be  produced  and  directed  at  wild- 
life in  their  natural  habitat  in  order  to  study  their  startle  response. 

In  order  to  illustrate  the  type  of  wave  that  must  be  produced 
by  a sonic-boom  simulator,  and  also  for  future  reference,  an  idealized 
overpressure  signature  of  a sonic  boom  is  sketched  in  Fig.  1.  The 
more  important  parameters  which  are  commonly  used  to  describe  the  various 
parts  of  the  signature  include  the  peak  overpressure,  rise  time,  dura- 
tion (or  wave  length),  and  wave  form  which  may  vary  somewhat  from  the 
ideal  'N'  shape.  It  is  worth  noting  that  respective  values  of  peak  over- 
pressure, duration  and  rise  time  are  100  N/m?,  300  ms  and  1 ms  for  a 
typical  sonic  boom  from  a current  SST  aircraft  and  also  from  large  mili- 
tary bomber  supersonic  aircraft.  In  the  case  of  shorter  military  fighter 
supersonic  aircraft  only  the  duration  is  significantly  different,  being 
correspondingly  shorter  at  about  100  ms. 

The  analytical  and  experimental  work  given  in  this  report  on 
the  wave  motion  in  low-pressure-ratio  rectangular  and  pyramidal  shock 
tubes  utilizing  different  driver  and  channel  gases  is  a natural  contin- 
uation of  previous  sonic-boom- simulati on  work  at  UTIAS.  Also,  the 
present  work  is  a direct  extension  of  sonic-boom-simulation  research  on 
shock  tubes  made  in  England,  France  and  the  United  States.  These  state- 
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merits  are  elaborated  on  in  the  following  part  of  this  introduction. 

A number  of  geometrically  different  shock  tubes  have  been  constr- 
ucted and  tested  by  English,  French,  American  and  Canadian  researchers  to 
assess  their  capability  of  producing  a good  'N-wave'  to  simulate  the  sonic 
boom  (Fig.  1).  Some  of  the  more  important  shock  tubes  are  illustrated 
schematically  in  Fig.  2,  and  the  developmental  and  research  efforts  made  by 
various  researchers  on  such  shock  tubes  are  summarized  in  Table  1.  Each 
shock  tube  shown  in  Fig.  2 has  been  named  according  to  the  geometrical  shape 
of  its  driver  and  channel.  The  adjectives  'rectangular'  and  'pyramidal' 
arise  quite  naturally  owing  to  the  geometry  of  the  driver  and  channel  of 
shock  tubes  commonly  used  in  the  laboratory.  More  generally,  however,  in 
this  report  rectangular  will  refer  to  any  duct  where  the  cross-sectional 
area  is  constant  with  distance  (e.g.,  cylindrical),  and  pyramidal  will 
refer  to  any  duct  where  the  cross-sectional  area  increases  directly  with  the 
square  of  the  distance  from  the  center  of  symmetry  (e.g.,  conical). 

For  shock  tubes  commonly  used  for  sonic-boom  simulation  purposes 
the  driver  and  channel  gases  are  both  normally  air  at  atmospheric  tempera- 
ture. The  driver  air,  however,  is  initially  at  a slightly  higher  pressure 
than  the  atmospheric  pressure  air  in  the  channel.  Consequently  the  pres- 
sure ratio  for  the  driver  and  channel  gases  separated  initially  by  a thin 
diaphragm  is  very  nearly  equal  to  unity.  On  breaking  the  diaphragm  in 
such  a low-pressure-ratio  shock  tube  the  ensuing  wave  motion  results  in  a 
very  weak  shock  wave  or  simulated  sonic  boom  moving  in  the  channel  gas. 

Closed-form  solutions  for  the  wave  motion  in  low-pressure-ratio 
rectangular , pyramidal,  pyramidal- rectangular , pyramidal-pyramidal  and 
rectangular-pyramidal  shock  tubes  which  use  identical  driver  and  channel 
gases  have  been  obtained  previously  (Refs.  5 and  6).  The  rectangular  shock 
tube  produces  a constant- amplitude  pulse  as  illustrated  by  the  overpres- 
sure signature  sketched  in  Fig.  3a-  In  the  case  of  the  pyramidal  shock  tube 
an  N-shaped  pulse  is  produced  as  shown  in  Fig.  3h.  The  other  pyramidal- 
rectangular,  pyramidal-pyramidal  and  rectangular-pyramidal  shock  tubes  each 
produces  a distorted  N-wave  followed  by  additional  waves  or  disturbances 
as  illustrated  in  Fig.  3c,  3d  and  3e  respectively.  Obviously  the  pyramidal 
shock  tube  is  best  for  producing  a good  N-wave  for  sonic-boom  simulation 
purposes.  Note  that  the  rapid  pressure  rise  across  both  the  front  and  rear 
shocks  of  the  N-wave  equals  ApQr0/2r  and  the  duration  of  the  N-wave  is  2r0/a, 
where  the  respective  symbols  Ap0,  a,  r and  rD  denote  the  pressure  difference 
across  the  diaphragm,  sound  speed  of  the  driver  gas,  radial  distance  measured 
from  the  driver  apex,  and  the  diaphragm  location  or  driver  length. 

A large  pyramidal  shock  tube  is  required  if  the  N-wave  produced 
in  the  channel  is  to  simulate  a full-scale  sonic  boom  from  an  SST  aircraft. 
Since  the  duration  of  the  N-wave  in  the  shock  tube  equals  2r0/a  and  the 
duration  of  a simulated  sonic  boom  has  to  be  approximately  300  ms,  the  cor- 
respDniing  length  of  the  shock-tube  (rQ)  needs  to  be  about  50  m.  Owing  to 
inadequate  breakage  of  a large  diaphragm  in  a low-pressure- ratio  shock  tube 
(Ref.  12),  the  largest  diaphragm  having  sufficiently  good  breaking  charac- 
teristics appears  to  be  about  1 nr  in  area.  The  maximum  divergence  angle  of 
a pyramidal  shock  tube  is  therefore  about  one-fiftieth  of  a radian.  Now  the 
interior  test  section  of  the  pyramidal  channel  must  be  at  least  3 m on  each 


2 


side  to  adequately  accommodate  a human  subject  without  unduly  blocking  the 
path  of  the  simulated  sonic  boom.  For  such  a test  section  size  and  diver- 
gence angle  the  resulting  length  of  the  pyramidal  shock  tube  is  an  incredible 
150  m.  In  practice  the  shock  tube  is  even  longer  because  a reflection 
eliminator  is  required  to  cover  the  open  end  of  the  channel.  For  the  sake  of 
interest  the  large  Franco-German  pyramidal  shock  tube  (Ref.  l4),  which  is 
189  m in  length,  is  depicted  in  Fig.  *4.  Note  that  different  diaphragm  sta- 
tions have  been  provided  such  that  the  N-wave  duration  can  be  varied  conveniently 
from  100  to  300  ms . 

Owing  to  the  large  size  required  of  a pyramidal  shock  tube  if  a 
full-scale  sonic  boom  is  to  be  simulated,  the  initial  construction  can  be 
a costly  endeavour.  If  a pyramidal-pyramidal  shock  tube  instead  of  a pyramidal 
shock  tube  could  be  used,  then  the  shock -tube  length  and  thus  the  initial  con- 
struction cost  might  be  reduced.  The  idea  behind  this  scheme  is  brought  out 
by  the  diagrams  given  in  Fig.  5.  In  order  to  maintain  the  same  sized  driver, 
diaphragm  and  test  section  of  a pyramidal  shock  tube  (Fig.  5a),  the  divergence 
angle  of  the  channel  is  increased,  thereby  producing  a shorter  shock  tube 
(Fig.  5b).  This  resulting  pyramidal-pyramidal  shock  tube,  however,  produces 
a distorted  N-wave  (Fig.  3d),  which  is  normally  undesirable  for  sonic-boom 
simulation  purposes.  In  fact  it  is  worth  mentioning  that,  based  on  the  ex- 
perimental results  obtained  from  prototype  pyramidal-pyramidal  shock  tubes 
(Ref.  12),  French  and  German  researchers  deemed  the  distortion  in  the  N-wave 
as  unacceptable,  and  they  consequently  constructed  the  large  Franco-German 
pyramidal  shock  tube  mentioned  previously  and  shown  in  Fig.  4. 

Another  possible  method  of  achieving  a shorter  and  thus  less 
costly  pyramidal  shock  tube  for  the  simulation  of  a full-scale  sonic  boom 
is  to  use  a driver  gas  which  has  a low  sound  speed.  (Although  the  channel 
gas  could  be  the  same  as  the  driver  gas,  it  would  probably  be  more  conven- 
ient to  use  air  in  the  channel.)  Because  the  predicted  duration  of  the 
N-wave  produced  in  the  shock  tube  is  2rQ/a,  if  a gas  having  a lower  sound 
speed  (a)  than  air  is  used  in  the  driver,  the  correct  full-scale  duration  of 
a sonic  boom  can  be  accomplished  with  a shorter- than-normal  driver  (length 
rQ)  and  hence  a shorter  pyramidal  shock  tube.  For  example,  the  use  of  an 
economical  gas  like  carbon  dioxide  (CO2)  instead  of  air  in  the  driver  and 
air  in  the  channel  would  result  in  a fairly  significant  reduction  of  22% 
in  shock-tube  length.  The  use  of  more  exotic  arid  expensive  dichlorodifluoro- 
methane  (C  CI2  F2)  > sulfur  hexafluoride  (SFg)  or  octofluorocyclobutane  (CI4 
Fg)  would  result  in  an  even  more  significant  reduction  in  shock- tube  length 
by  56%,  60%  or  67%  respectively. 

Another  possible  advantage  of  using  a driver  gas  having  a low 
sound  speed  is  to  improve  the  capability  of  an  existing  low-pressure-ratio 
pyramidal  shock  tube.  Most  of  the  existing  shock  tubes  have  been  designed 
to  be  short  to  minimize  construction  costs,  and  they  consequently  produce 
a short- duration  simulated  sonic  boom.  (Such  short-duration  booms  can  be 
used  for  studies  of  wave  diffraction  over  and  into  model  buildings,  wave 
propagation  over  reduced- scale  land  topologies,  and  certain  animal  response 
tests) . Since  atmospheric  temperature  air  is  normally  used  in  both  the 
driver  and  channel  the  standard  technique  of  changing  the  N-wave  duration 
for  new  tests  is  to  alter  the  driver  length  by  changing  the  diaphragm  loc- 
ation. As  existing  shock  tubes  are  normally  short  (less  than  15  m long)  an 
upper  limit  exists  on  the  driver  length  and  therefore  on  the  N-wave  duration. 
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Consequently,  the  only  means  of  further  increasing  the  N-wave  duration  and 
thus  extending  the  capability  of  existing  shock  tubes  is  to  use  a driver 
gas  which  has  a sufficiently  low  sound  speed.  For  example,  the  N-wave  dur- 
ation can  be  increased  markedly  by  a factor  of  1.3,  2.3,  2.5  or  3*0  when 
the  driver  gas  is  CO2,  CC12F2,  SF^  or  C4F8,  respectively,  instead  of  air. 

It  should  be  noted  that  the  N-wave  duration  can  be  increased  further  by 
lowering  the  temperature  of  the  driver  gas.  However,  this  method  is  not 
only  impracticable  but  the  increase  in  duration  is  not  very  significant, 
being  only  20$>  for  a decrease  in  driver-gas  temperature  from  300  to  200  K. 

The  technique  of  using  a driver  gas  with  a low  sound  speed  in  a 
pyramidal  shock  tube  to  increase  the  N-wave  duration  both  to  extend  the 
capability  of  an  existing  shock  tube  and  to  investigate  the  feasibility 
of  using  different  driver  and  channel  gases  in  a larger  sonic-boom  simula- 
tor was  originally  tested  in  England  (Ref.  9)*  Carbon  dioxide  was  used 
in  the  driver  of  a small  conical  shock  tube  (driver  length  of  60  cm, 
channel  length  of  lUO  cm)  and  showed  that  a good  N-wave  could  be  produced 
in  the  channel  air.  The  N-wave  duration  was  an  expected  1.3  times  longer 
than  that  for  the  case  when  air  was  used  in  both  the  driver  and  channel. 

To  obtain  an  even  longer  duration  N-wave  with  the  same  shock  tube  the 
use  of  a special  gas  having  the  trade  name  'Arcton'  was  proposed,  as  its 
low  sound  speed  is  only  kjf0  of  that  for  air.  Test  results  for  the  case 
of  Arcton  as  the  driver  gas  were  unfortunately  not  reported. 

The  British  researchers  (Ref.  9)  did  not  notice  any  significant 
distortion  in  the  overpressure  signature  of  the  N-wave  produced  in 
their  conical  shock  tube  when  carbon  dioxide  was  used  in  the  driver  and 
the  channel  gas  was  air.  However  experiments  at  UTIAS  showed  that  if  the 
sound  speeds  of  the  driver  and  channel  gases  were  radically  different 
(e.g.,  CCI2F2,  SFg  or  C4F8  being  the  driver  gas  with  air  in  the  channel) 
then  a noticeably  distorted  N-wave  would  be  produced.  This  distortion  in 
the  N-wave  is  similar  to  that  in  the  N-wave  produced  by  a pyramidal -pyra- 
midal shock  tube  which  used  identical  driver  and  channel  gases  (Fig.  3d) . 

If  the  undesirable  distortion  in  the  N-wave  is  significant  then  the  con- 
cept of  using  different  driver  and  channel  gases  in  a pyramidal  shock 
tube  either  to  increase  the  duration  of  an  N-wave  produced  in  an  existing 
short  shock  tube  or  to  enable  a short  and  inexpensive  shock  tube  to  pro- 
duce a simulated  full-scale  sonic  boom  is  not  valid.  This  will  be  verified 
subsequently  in  this  report  by  analytical  and  experimental  data. 

The  acoustic  analysis  for  the  wave  motion  in  a pyramidal  shock 
tube  utilizing  different  driver  and  channel  gases  was  found  to  be  quite 
complex  analytically.  Prior  to  obtaining  the  acoustic  solution  for  the 
pyramidal  shock  tube  a similar  but  simpler  acoustic  solution  was  obtained 
for  the  case  of  a rectangular  shock  tube  utilizing  different  driver  and 
channel  gases.  For  interest  and  completeness  both  of  these  acoustic 
analyses  are  given.  The  simpler  analysis  for  the  rectangular  shock  tube 
is  presented  first  in  Chapter  2 and  the  more  complex  analysis  for  the 
pyramidal  shock  tube  follows  in  Chapter  3*  Experimental  data  for  both 
the  rectangular  and  pyramidal  shock  tubes  are  also  given,  verifying  the 
analyses. 

It  is  worth  mentioning  that  the  respective  acoustic  analyses  for 
the  rectangular  and  pyramidal  shock  tubes  are  applicable  directly  to  weak 
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planar  and  spherical  explosions  of  finite  size.  Consequently  the  acoustic 
analyses  are  of  fundamental  importance  vo  the  understanding  of  the  wave 
motion  of  planar  and  spherical  explosions,  for  which  the  explosion  gas 
differs  from  the  ambient  gas . 


2.  RECTANGULAR  SHOCK  TUBE 

The  rectangular  shock  tube,  which  is  depicted  in  Fig.  2a,  con- 
sists essentially  of  a constant  cross-sectional  area  driver  and  channel 
which  are  joined  at  the  diaphragm  station.  The  diaphragm  initially  sep- 
arates a normally  atmospheric  pressure  channel  gas  from  a slightly  higher 
pressure  driver  gas.  For  generality  of  the  analysis  let  the  initially 
quiescent  driver  and  channel  gases  be  not  only  different  but  also  have  dif- 
ferent temperatures . Before  the  diaphragm  is  broken  the  appropriate  initial 
conditions  of  the  driver  and  channel  gases  for  the  analysis  can  therefore 
be  summarized  mathematically  as  follows. 


Channel  (x  > x ) : 
' o 

Ap  = 0 

(2.1) 

o 

II 

(2.2) 

Driver  ( 0 < x < x ) : 

o 

o 

3 

II 

OJ 

3 

(2.3) 

Ai^  = 0 

(2.4) 

The  respective  symbols  Ap,  Au,  and  Apc  denote  overpressure  with  respect  to 
atmospheric  pressure,  perturbation  flow  velocity  or  particle  velocity,  and 
pressure  difference  across  the  diaphragm,  whereas  the  respective  subscripts 
1 and  2 refer  the  appropriate  symbols  to  the  channel  and  driver  gases . The 
distance  x is  measured  along  the  shock  tube  starting  from  the  closed  end  of 
of  the  driver  as  shown  in  Fig.  2a,  and  the  diaphragm  is  located  at  distance 
xQ,  thereby  making  the  driver  length  equal  to  xQ. 

For  the  preceding  initial  conditions  the  ensuing  wave  motion  in 
the  rectangular  shock  tube  after  the  diaphragm  is  broken  is  rather  complex. 
This  wave  motion  can  be  depicted  conveniently  on  a time-distance  diagram, 
as  shown  in  Fig.  6,  where  the  locus  of  each  wave  front  has  been  drawn.  On 
breaking  the  di  aphragm  in  the  shock  tube  the  resulting  rapidly  expanding 
driver  gas  produces  a weak  shock  wave  (g^)  in  the  channel  gas.  Simultan- 
eously a weak  rarefraction  wave  (f^)  moves  into  the  driver  gas  and 
eventually  reaches  and  reflects  from  the  closed  end.  This  reflected  wave 
(h^)  propagates  through  the  driver  gas  and  encounters  the  contact  surface 
or  interface  of  the  driver  and  channel  gases.  Because  the  contact- surface 
displacement  is  negligible  in  comparison  with  the  driver  length,  the  con- 
tact- surf ace  path  is  shown  in  Fig.  6 singly  as  a vertical  dashed  line. 

Owing  to  the  change  in  specific  impedance  (product  of  the  density  and  sound 
speed)  across  the  contact  surface  the  reflected  wave  (h^)  is  partially 
transmitted  to  the  channel  (gj>)  and  partially  reflected  back  into  the  driver 
(f2).  The  subsequent  oscillatory  wave  motion  between  the  contact  surface 
and  closed  end  of  the  driver  creates  a sequence  of  waves  in  the  channel. 

The  integrated  result  of  all  of  these  individual  waves  (gi,  i = 1,  2,  ..., 
oo)  yields  the  total  wave  in  the  channel. 
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For  a low-pressure- ratio  rectangular  shock  tube  for  which  the 
pressure  difference  across  the  diaphragm  Apo  is  much  less  than  the  absolute 
pressure  of  the  channel  gas,  the  individual  and  total  waves  in  the  driver 
and  channel  are  relatively  weak.  The  wave  motion  can  therefore  be  described 
adequately  by  using  well-known  acoustic  theory.  The  one-dimensional  planar 
wave  equations  which  govern  the  acoustic  wave  motion  in  the  channel  and 
driver  gases  are  given  below. 


Channel  (x  > xo): 


Driver  (0  < x < xq) : 


S2o1 

2 

S2o1 

dt2 

= *1 

Sx2 

82t>2 

2 

82<t>2 

dt2 

= *2 

Sx2 

(2.5) 


(2.6) 


The  respective  symbols  , a and  t denote  total  velocity  potential,  sound 
speed  and  time. 


For  the  analysis  it  is  convenient  to  express  the  total  velocity 
potential  of  the  channel  gas  (<t>i)  and  the  driver  gas  (<>2)  in  a more  basic 
form  consisting  of  the  sum  of  velocity  potentials  of  individual  waves  in 
the  channel  and  driver  respectively.  The  new  forms  of  <tq  and  <t>2 , and 
related  expressions  for  overpressure  and  particle  velocity,  are  summarized 
below. 


Channel  (x  > xq): 


*1  =X[gi(l)  ” 2(i"l)al/a2  + 1} 

i=l 

Api  = -pi  w 

M1  Apx 

Aui = sir = z[p7 

Driver  (0  < x < xq): 

Ap  t °* 

+ £[fi(l)  Hfo-2i+l)  +^0)  H{p-2i+})j 
2 i=l 

d<>2 

Ap2  = ’P2  W 

A 

Au2  = sr 


(2.7) 

(2.8) 

(2.9) 


(2.10) 

(2.11) 

(2.12) 
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In  the  above  expressions:  p denotes  density;  the  respective  nondimen- 

sional  variables  I,  i},  P equal  (axt  - x)/xQ,  (agt  + x)/xQ  and  (a2t  - x)/xQ; 
Si(0»  fi(f|)  and  hi(P)  denote  the  velocity  potentials  of  the  itb  distur- 
bances as  illustrated  in  Fig.  6;  and  H(Tj.}  denotes  the  unit  step  function 
which  equals  zero  prior  to  the  arrival  of  the  ith  wave  (Ti  < 0)  and  equals 
unity  after  the  passage  of  the  ith  wave  front  (T*  > 0).  The  form  of  the 
total  velocity  potentials  0j.  (Eq.  2.7)  and  02  (Eq.  2.10)  are  quite 
general  but  not  arbitrary,  as  they  have  been  chosen  to  satisfy  both  their 
respective  wave  equations  (Eqs.  2.5  and  2.6)  and  their  associated  initial 
conditions  (Eqs.  2.1  and  2.2  and  Eqs.  2.3  and  2.4).  As  an  example  of 
satisfying  the  initial  conditions  the  term  -^0t/p2  in  Eq.  2.10  accounts 
for  the  pressure  difference  across  the  diaphragm  (£f>0)  for  time  t less  than 
zero. 


The  presently  unknown  velocity  potentials  gi(l),  fi(ij)  and  h^(P) 
can  now  be  determined  by  evaluating  the  effects  of  two  different  boundaries 
on  the  wave  motion.  The  first  boundary  to  be  considered  is  the  closed 
end  of  the  driver.  At  this  stationary  boundary  the  particle  velocity  of 
the  driver  gas  must  be  zero  for  all  time.  By  setting  the  particle  velocity 
of  the  driver  gas  AU2  (Eqs.  2.10  and  2.12)  equal  to  zero  at  this  boundary 
where  distance  x equals  zero  the  following  intermediate  result  can  be 
obtained. 


^ [^q(ij)  H(T)-2i+l]  - h^(P)  H(P-2i+l) 
i=l  L 


= 0 


(2.13) 


The  prime  (')  denotes  differentiation  of  the  variable  with  respect  to  the 
argument  given  in  the  following  brackets.  Owing  to  the  mathematical  struc- 
ture of  this  expression  (Eq.  2.13)  it  will  in  general  be  identically  zero 
for  all  time  if  and  only  if  certain  terms  having  equivalent  step  functions 
cancel  exactly.  As  t|  and  p both  equal  a2t/x0  at  the  boundary  x equals  zero, 
the  lengthy  expression  of  Eq.  2.13  can  consequently  be  expressed  in  the 
following  equivalent  and  more  convenient  form. 


fcj(P)  = ^U)  i = 1,  2, 


(2.14) 


On  integrating  each  one  of  these  first-order  differential  equations  with 
respect  to  time,  and  after  setting  the  constants  of  integration  equal  to 
zero  as  they  are  arbitrary,  the  following  desired  result  can  be  obtained. 


h^P)  ■ ^(il)  i = 1,  2,  ...,  oo 


(2.15) 


These  results  illustrate  that  the  particle  velocity  of  the  ifch  wave  having 
a velocity  potential  fi(rj)  is  countered  exactly  at  the  closed  end  of  the 
driver  by  the  particle  velocity  of  the  corresponding  ith  reflected  wave 
having  a velocity  potential  hj^P),  thereby  having  a net  particle  velocity 
of  zero  at  this  stationary  boundary. 
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The  second  boundary  to  be  considered  is  the  interface  of  the  driver 
and  channel  gases.  At  this  contact  surface  the  overpressure  and  particle 
velocity  are  taken  to  be  continuous,  and  the  effects  of  diffusion,  heat 
transfer  and  turbulent  mixing  are  neglected.  The  matching  of  the  driver  and 
channel  gas  overpressures  and  particle  velocities  at  the  contact  surface  can 
be  difficult  mathematically  because  the  contact  surface  is  not  stationary 
nor  is  its  motion  known  a priori.  To  circumvent  this  difficulty  the  effects 
of  the  contact- surface  motion  are  assumed  to  be  negligible  and  the  matching 
procedure  is  performed  at  the  diaphragm  station  - a fixed  or  stationary 
location  xQ. 


For  the  first  step  in  the  matching  procedure  at  the  diaphragm 
station  (x  = xG)  the  overpressures  of  the  channel  gas  £pq  (Eqs.  2.7  and 
2.8)  and  the  driver  gas  Ap2  (Eqs.  2.10  and  2.11)  are  equated,  yielding 
the  following  intermediate  result. 


= -fa  X 
^O  o 


y>l)Tk<0  H{!  - 2(i-l)a1/a2  + 1} 
i=l 

00 

,^[q(n)  H(ri-2i+l)  + h’(p)  H{P-2i+l)  (2.16) 


+ a2P2 


For  this  boundary  (x  = xQ)  £ equals  (aqt  - x0)/x0,  p equals  (a2t  - x0)/x0 
and  i)  equals  (a2t  + x0)/x0.  Certain  terms  in  Eq.  2.16  consequently  have 
equivalent  step  functions,  which  can  be  grouped  accordingly,  such  that 
Eq.  2.16  can  be  rewritten  in  the  following  equivalent  and  more  convenient 
form. 


alpl  gl(5)  = -^oxo  + a2p2  1 = 1 (2.17) 

alpl  gi(^  = a2p2  hi-l^^  + a2P2  fi(T|^  i = 2,  3,  (2.18) 


The  first  expression  (Eq.  2.17)  shows  that  at  the  contact  surface  the  over- 
pressure of  the  first  (shock)  wave  in  the  channel  (gq)  matches  exactly  the 
sum  of  the  initial  overpressure  in  the  driver  (£p0)  and  the  overpressure  of 
the  first  (rarefaction)  wave  in  the  driver  (fq).  The  second  expression 
(Eq.  2.18)  shows  that,  for  the  i^h  wave  interaction  at  the  contact  surface 
(see  Fig.  6),  the  sum  of  the  overpressures  of  the  incident  wave  (hq_q)  and 
reflected  wave  (fq)  matches  exactly  the  overpressure  of  the  transmitted  wave 
(gi)* 

The  second  and  final  step  of  the  matching  procedure  at  the  diaphragm 
station  (x  = xQ)  is  to  equate  the  particle  velocity  of  the  channel  gas  Auq 
(Eqs.  2.7  and  2.9)  to  that  of  the  driver  gas  Au2  (Eqs.  2.10  and  2.12),  giving 
the  following  intermediate  result. 
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* 


00^ 

]T  [gi(0  HU  - 2(i-l)a1/a2  + 1} 

L=1 

00 

£ [h'(p)  H{p-2i+l)  - f'(T])  H(n-2i+l} 


(2.19) 


Like  the  previous  result  for  the  overpressure  (Eq.  2.16)  certain  terns  in  Eq. 
2.19  have  equivalent  step  functions,  which  can  be  grouped  accordingly,  and 
Eq.  2.19  rewritten  as  shown  below. 


g£  (!)  --fi(l)  i=l  (2.20) 

g'  (i)  = h'^U)  - q(t|)  i=2,  3,  ...,  oo  (2.21) 

The  first  expression  (Eq.  2.20)  shows  that  at  the  contact  surface  the  par- 
ticle velocity  of  the  first  (shock)  wave  in  the  channel  (g^)  matches  exactly 
that  of  the  first  ( rarefaction)  wave  in  the  driver  ( f]_) . The  second  ex- 
pression (Eq.  2.21)  shows  that,  for  the  i^*1  wave  interaction  at  the  contact 
surface  (see  Fig.  6),  the  sum  of  the  particle  velocities  of  the  incident 
wave  (H^_^)  and  reflected  wave  (f^)  matches  exactly  the  particle  velocity  of 
the  transmitted  wave  (gi). 

There  is  actually  a third  boundary  which  is  the  end  of  the  shock- 
tube  channel.  However,  this  boundary  need  not  be  considered  here  as  it  has 
already  been  illicitly  assumed  that  the  channel  is  either  infinitely  long 
or  terminated  by  a perfect  reflection  eliminator  such  that  no  reflected  or 
other  waves  arise  from  this  boundary. 

Now  that  the  matching  procedure  has  been  con$?leted  the  resulting 
equations  can  be  used  to  obtain  final  expressions  for  the  i^h  velocity 
potentials  g^U),  fid)  and  hiO)  . From  Eqs.  2.17  and  2.20,  which  are 
simultaneous  expressions  in  g£(|)  and  f£d) , one  can  obtain  g{ ( I)  and  fl(t)) 
explicitly.  In  a similar  manner,  from  Eqs.  2.18  and  2.20,  g |(i)  and  f{(t)) 
for  i greater  than  unity  can  be  obtained  in  terms  of  h^_^(3).  These  results 
are  summarized  below. 


g{U) 


gj(!) 

. 

q<i> 


-'Vo 

i=l 

Vi  + Vs 

(2.22) 

Vi  * Vs  hl-1<e) 

i=2,  3, 

• • • j 00 

Ap  x 
o o 

i=l 

Vi  * Vs 

(2.23) 

- ^ I V2  K l(s) 

i=2,  3, 

00 

• • • > 

Vi  + Vs  i‘1 
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Then,  in  a straight-forward  procedure  using  Eqs.  2.14,  2.22  and  2.23, 
explicit  expressions  can  be  obtained  for  g!(i),  f 7 ( t])  and  h'  (P).  These 
expressions  can  be  easily  integrated  to  yield  the1final  results  for 
gi(0>  and  h^P),  which  are  given  below. 


g±(0 

^(0 


^(3) 


-Ap  x ? 
^o  o 


- vfr 


a2p2 


a1P1  + a2P2 

aj.Pj.  ~ *2°2 

Vl  + a2P2 


(Vl  + *2°2' 

l-1  *0  1 

alPl  + a2P2 


Vl  “ a2P2  l1"1  Apo  Xq  P 

[ a1p1  + a2P2  a^P1  + a2P2 


i=l 


Vl  - a2p2  I1'2  2a2p2  ^o  Xo  5 , 

■JT  i=2>  3i 


(2.24) 


i=l,  2, . . .,c 


i=l,  2,...,c 


(2.25) 

(2.26) 


Note  that  the  constants  of  integration  have  all  been  set  equal  to  zero 
because  they  are  arbitrary,  and  g’(l),  f ’ ( tj)  and  h.’ (P)  can  be  easily  recov- 
ered by  differentiation.  Also,  these  final  expressions  for  g.(|),  f.(r)) 
and  h.  (P)  are  valid  away  from  the  two  boundaries  x equal  to  zero  and^x 
as  x In  t,  T)  and  P need  no  longer  be  restricted  to  zero  or  x 
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The  solution  for  the  wave  motion  in  a rectangular  shock  tube  util- 
izing different  driver  and  channel  gases  has  now  been  obtained.  The  final 
results  for  the  ith  velocity  potentials  g .(£),  f . (tj)  and  h (p)  (Eqs.  2.24, 
2.25  and  2.26)  can  be  substituted  into  the  expressions  forthe  total  velocity 
potentials  <t>,  and  <t>2  (Eqs.  2.7  and  2.10).  Then  the  two  respective  expres- 
sions for  the  overpressure  and  particle  velocity  of  the  channel  gas  follow 
from  Eqs.  2.8  and  2.9,  and  those  for  the  driver  gas  follow  from  Eqs.  2.11  and 
2.12. 


It  can  be  seen  from  the  solution  of  the  wave  motion  in  the  rectan- 
gular shock  tube  that  the  specific  impedances  of  the  channel  gas  (a-p..) 
and  the  driver  gas  (a2P2)  play  an  important  role.  In  fact,  the  important 
parameter  is  the  ratio  of  these  specific  impedances.  This  ratio  is  denoted 
by  the  symbol  * and  it  can  be  expressed  in  the  following  alternate  forms 
by  using  the  equation  of  state  (p=pRT)  and  the  isentropic  sound- speed  rela- 
tion (a2  = 7p/p  « 7r  T/M) . 


Vl 

7ia2 


r Tr 


^ Tll* 


l >1 


(2.27) 


The  respective  symbols  y,  p,  R,  T,  RQ  and  M denote  the  specific  heat  ratio, 
pressure,  gas  constant,  temperature , uni versed  gas  constant  (R0  = MR)  and 
molecular  weight.  It  is  worth  mentioning  that  the  impedance  ratio  9 equals 
unity  not  only  for  the  simplest  case  of  identiced  driver  and  channel  gases 
having  the  same  temperatures  but  also  for  the  specific  case  of  different 
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driver  and  channel  gases  provided  that  y VL/t  §quals  y^  NU/t^,  . Further- 
more, the  impedance  ratio  b equals  ( y 2 MT)2  for  the  very  important 

case  of  different  driver  and  channel  gases  having  equal  temperatures,  and 
b equals  (T^/T  )?  for  the  not  so  practicable  case  of  identical  driver  and 
channel  gases  having  different  temperatures. 

A set  of  values  of  the  specific  impedance  ratio  b (agP2/aiPi  = 
(72  Mo/7o  if  T^=T2)  was  compiled  for  certain  interesting  different 

combinations  of  driver  and  channel  gases  having  identical  temperatures, 
and  they  are  given  in  Table  2.  The  smallest  values  of  b occur  when  a 
light  driver  gas  such  as  hydrogen  (H2)  or  helium  (He)  is  used  in  conjunc- 
tion with  a heavy  channel  gas  such  as  dichlorodifluoromethane  (CCI2  F2) , 
sulfer  hexafluoride  (SF6)  or  octafluorocyclobutane  (C4  Fg) . On  the  other 
hand  the  largest  values  of  b occur  when  the  driver  gas  has  a high  molec- 
ular weight  and  the  channel  gas  has  a low  molecular  weight.  For  nearly 
equivalent  molecular  weights  of  the  driver  and  channel  gases  the  specific 
inpedance  ratio  is  always  nearly  equal  to  unity. 


The  wave  that  propagates  in  the  channel  gas  is  one  of  the  most 
important  features  of  the  wave  motion.  For  this  reason  the  signature  of 
pressure  signature,  as  derived  from  Eqs.  2.7,  2.8,  2.24  and  2.27,  is 
given  below: 


APi-  §E|  Hts*u 

i-2 


i-2  a 

H{|  - 2(i-l)  -±  + 

72 


1) 


(2.28) 


After  the  arrival  of  the  first  shock  wave  (gi)  at  time  (x-Xo)/a^  at  a fixed 
distance  x in  the  channel,  each  successive  wave  (g.^)  arrives  at  time  (x-xQ)/ 
a^  + 2(i-l)x0/a£i,  or  after  equal  time  intervals  of  2x0/ag.  Each  successive 
wave  has  an  initial  overpressure  of  zero  and  features  a sudden  change  in 
overpressure  at  its  wave  front  to  a constant  value  thereafter.  The  first 
wave  has  an  amplitude  of  Ap0/{l+s) , and  each  of  the  following  waves  has  a 
successively  smaller  amplitude  by  virtue  of  the  factor  (1-b)/(1+b)  raised  to 
the  power  i-2.  Because  each  successive  wave  is  superposed  on  the  previous 
ones  the  total  overpressure  signature  consists  of  an  infinite  sequence 
of  constant  overpressure  segments,  each  segment  having  a duration  of  2x.q/bq. 
The  sum  of  terms  giving  each  segment  of  the  signature  has  the  form  of  a 
finite  geometric  progression.  Consequently  the  following  geometric  sequence 
can  be  derived  for  the  amplitude  of  the  i^h  segment. 


li-1 


Apo 

1 + B 


i=l, 


(2.29) 


This  expression  is  more  convenient  to  use  when  constructing  the  overpressure 
signature  than  Eq.  2.28. 

Five  different  predicted  overpressure  signatures  of  the  wave  moving 
in  the  channel  are  shown  in  Fig.  7,  corresponding  to  specific  inpedance  ratios 
(a)  of  0.2,  0.5,  1.0,  2.0  and  5*0.  Although  these  irpedance  ratios  represent 
general  cases  of  different  driver  and  channel  gases  having  different  tempera- 
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tures,  various  important  combinations  of  equivalent-temperature  driver  and 
channel  gases  which  have  such  impedance  ratios  and  would  thus  produce  such 
signatures  can  be  readily  deduced  from  the  information  given  in  Table  2. 

The  top  two  signatures  in  Fig.  7 for  the  case  of  a light  driver  gas  and  a 
heavy  channel  gas  (low  « values  of  0.2  and  0.5)  exhibit  a series  of  very 
noticeable  descending  ’steps'  and  have  relatively  high  peak  overpressures 
(0.83Ap  and  0.67Apo)  for  a given  pressure  difference  across  the  dia- 
phragm (Ap0) . In  general  the  duration  of  each  constant  overpressure 
segment  (2x0/a2)  would  be  relatively  short  because  the  sound  speed  of  a 
light  or  heated  driver  gas  (ag)  in  relatively  high  (see  the  information 
in  Table  2 on  the  sound  speed  of  different  driver  gases) . The  center 
signature  for  the  case  when  * has  the  value  of  unity  exhibits  a single  step 
having  an  amplitude  of  0.5Apo  and  a duration  of  2xQ/a2 . This  signature 
corresponds  to  the  well-known  case  of  identical  driver  and  channel  gases 
having  the  same  temperature.  However,  it  can  also  correspond  to  the 
specific  case  of  different  driver  and  channel  gases  provided  72  %Ai 
equals  72  1^/T2-  The  bottom  two  signatures  for  the  case  of  a heavy  driver 
gas  and  a light  channel  gas  (high  values  of  b of  2 and  5)  exhibit  a series 
of  steps  that  alternate  markedly  in  amplitude  from  positive  to  negative 
values  as  they  diminish  in  absolute  magnitude.  The  relatively  low  peak 
overpressure  of  these  two  signatures  (0.33Apo  and  0.17£po)  in  comparison 
to  the  pressure  difference  across  the  diaphragm  (Ap0)  illustrates  the 
inefficiency  of  using  a heavy  driver  gas  to  produce  a relatively  strong  wave 
in  a light  channel  gas.  Note  that  these  two  signatures  will  have  rela- 
tively long  durations  because  the  sound  speed  of  the  heavy  driver  gases  will 
generally  be  low. 

The  particle-velocity  signature  of  the  wave  moving  in  the  channel 
gas  has  the  same  features  as  the  overpressure  profile,  because  the  particle 
velocity  of  this  planar  wave  is  directly  proportional  to  the  overpressure 
(Au-^  = Ap-^/a^Pl).  The  preceding  remarks  concerning  the  overpressure  profile 
therefore  apply  equally  well  to  the  particle-velocity  signature. 

The  motion  of  the  channel  gas  and  contact  surface,  owing  to  the 
wave  motion  in  the  rectangular  shock  tube,  can  be  determined  in  the  fol- 
lowing manner.  The  particle-velocity  signature  like  the  overpressure 
profile  consists  of  a sequence  of  descending  or  alternating  constant  over- 
pressure steps  or  segments.  The  particle  velocity  (Au^  - Zpi/a^pq) 
associated  with  the  i^*1  segment  of  its  profile  follows  from  Eq.  2.29  for  the 
overpressure  and  is  given  below. 

1__  ..  1-g  1 1-1  i = 1,  2,...,  00  (2.30) 

1 + a 1+h  a^^ 

During  each  constant  particle-velocity  segment  of  the  profile,  which  has  a 
fixed  duration  of  2x0/s^,  the  particle  path  will  be  a linear  function  of 
time.  The  displacement  of  a fluid  particle  by  the  i^*1  segment  of  the  wave 
is  thus  equal  to  the  result  of  Eq.  2.30  multiplied  by  2-Kq/bq.  The  accumu- 
lated displacement  of  a fluid  particle  after  n successive  segments  of 
duration  2x  /ag,  or  at  successive  times  2nxo/a2  measured  from  the  wave 
front,  can  be  determined  by  addition,  and  this  result  is  shown  below. 
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(2.31) 


1 - g 

1 + B 


-li-l 


00 


Because  this  sum  is  in  the  form  of  a geometric  progression,  Eq.  2.3.1  can 
he  easily  expressed  in  an  alternate  and  shorter  form.  By  noting  that 
ak  equals  7^p0/p2  or  72p  /pg  to  first  order,  the  final  expression  for 
E^.  2.31  takes  the  following  form. 


££cL*a 

72  Px 


n = 1,  2,, 


(2.32) 


In  summary,  this  expression  gives  the  displacement  of  a fluid  particle  x-x* 
from  its  initial  location  x*  at  successive  times  2nx0/a 2 measured  from  the 
wave  front.  The  displacement  during  each  time  interval  2x0/a^  is  a linear 
function  of  time. 

The  wave  moving  in  the  channel  gas  induces  the  same  displacemert 
x-x*  for  each  fluid  particle  because  its  predicted  signature  and  amplitude 
are  invariant  with  distance.  Consequently,  the  contact- surface  displacement 
x-x0  is  also  given  by  Eq.  2.32.  The  total  displacement  of  any  fluid 
particle  and  also  the  contact  surface  for  large  times  (n->»)  becomes  in- 
dependent of  the  specific  impedance  ratio  (a)  and  tends  to  a constant  value 
of  Ap  x0/72  p-^.  This  total  displacement  is  normally  small  when  compared 
with  the  driver  length  (xQ).  For  exanqple,  even  if  the  pressure  difference 
across  the  diaphragm  (Apc)  is  quite  high  at  p^/20  (one  twentieth  of  an 
atmosphere)  and  72  equals  1.1*  for  a diatomic  driver  gas,  the  total  dis- 
placement is  only  x0/28  or  of  the  driver  length.  Note  that  the 

original  assumption  for  the  analysis  that  the  contact- surface  displacement 
was  negligible  compared  with  the  driver  length  was  therefore  very  reasonable. 

Three  different  particle  paths  corresponding  to  values  of  the 
specific  impedance  ratio  ■ of  0.2,  1.0  and  5.0  are  shown  superposed  on 
the  time-distance  diagram  in  Fig.  8.  The  displacement  x-x*  of  a fluid 
particle  frcm  its  initial  location  x„  has  been  exaggerated  for  clarity. 

For  the  ctse  when  ■ equals  0.2  the  displacement  increases  monotomically  to 
its  maximum  value  of  ApQ  xq/72  P-.  • This  particle-path  behaviour  is  typical 
for  the  case  when  ■ lies  in  the  range  from  zero  to  unity,  or  when  a light 
driver  gas  is  used  in  conjunction  with  a heavy  channel  gas.  For  a driver 
and  channel  gas  having  a specific  impedance  ratio  of  unity  the  displace- 
ment increases  linearly  to  its  maximum  value  Ap  xq/72  p^  in  a time 
interval  of  2xo/a 2,  after  which  the  displacement  is  constant.  For  cases 
when  the  specific  'impedance  ratio  is  greater  than  unity  the  displacement 
alternately  increases  and  decreases  linearly  as  the  final  value  £pQ 
7p  Pt  is  approached,  as  illustrated  in  Fig.  8 by  the  particle  path  for 
which  b equals  5. 

The  final  displacement  of  the  contact  surface  x-xQ,  which  equals 
Ap  xj 72  p.,  can  be  derived  in  an  alternate  manner,  thereby  providing  a 
valuable  independent  check  on  the  results  of  the  preceding  analysis  for  the 
rectangular  shock  tube.  Let  the  driver  gas  have  an  initial  pressure  of 
p.  + Ap  and  an  initial  density  given  by  m /Ax  , where  m is  the  total  mass 
or  the  °driver  gas,  AxQ  is  the  driver  volume  and  A is  th?  cross-sectional 
area  of  the  shock  tube.  Now  let  the  driver  gas  expand  i sent Topically  (as 
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it  would  for  acoustic  wave  motion)  to  a final  density  of  mo/Ax  and  a 
final  pressure  of  which  is  the  initial  pressure  of  the  channel  gas. 
The  following  isentropic  expression  relating  the  initial  and  final 
pressures  and  densities  can  therefore  be  applied. 


Pi  + APq  ! 

r VAxoy 

2 _ 

r x + x - x -i 
| o o 

*d 

H 

1 

L J 

1 x 

L O J 

(2.33) 


Consequently,  to  first  order,  the  final  contact- surface  displacement 
x-x0  equals  Ap0  x0/?2  Pi,  in  exact  agreement  with  the  result  of  the 
previous  acoustic  analysis. 


A further  interesting  and  independent  check  on  the  acoustic 
analysis  for  the  rectangular  shock  tube  is  provided  by  conventional 
shock-tube  theory,  which  can  be  found  in  Ref.  17.  For  a constant- area 
shock  tube  utilizing  perfect  driver  and  channel  gases  having  a small  or  large 
pressure  difference  across  the  diaphragm  (ApQ),  the  following  conven- 
tional shock-tube  equation  can  be  used  to  predict  the  peak  overpressure 
(Ap. ) of  the  constant- amplitude  shock  wave  in  the  channel  gas. 

-27^ 
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From  this  equation  it  can  be  easily  shown  that,  when  dp  is  much  smaller 
that  p.,  the  peak  overpressure  of  the  first  shock  wave  (Api)  equals 
Ap0/(l+  b),  in  exact  agreement  with  the  corresponding  result  of  the 
previous  acoustic  analysis.  It  is  worth  mentioning  that  although  the 
conventional  shock-tube  equation  is  valid  for  large  as  well  as  small 
pressure  differences  across  the  diaphragm,  the  equation  is  limited  in  that 
it  can  predict  the  amplitude  of  only  the  first  but  most  inportant  part 
of  the  wave  in  the  channel  gas.  Admittedly,  additional  shock- tube  theory 
can  be  applied  oo  predict  the  complete  wave  motion  for  all  time  in  the 
constant-area  shock  tube.  However,  the  procedure  is  prohibitively 
difficult.  Consequently,  the  acoustic  solution  for  the  complete  wave 
motion  in  the  shock  tube,  even  though  it  is  restricted  to  acoustic  waves, 
is  very  valuable. 


Some  overpressure  measurements  were  made  in  ’the  channel  of  a 
constant-area  shock  tube,  in  order  to  provide  an  experimental  check 
on  the  validity  of  the  acoustic  analysis  for  the  rectangular  shock  tube. 
Three  of  the  measured  overpressure  signatures  are  presented  in  Fig.  9- 
For  the  profile  shown  in  Fig.  9a  the  driver  and  channel  gases  were  both 
air  having  the  same  temperature.  The  profile  shown  in  Fig.  9b  corresponds 
to  the  case  of  helium  in  the  driver  and  equivalent-tenperature  air  in  the 
channel.  For  the  opposite  case  of  air  in  the  driver  and  helium  in  the 
channel,  the  signature  appears  in  Fig.  9C*  Other  pertinent  data  are 
given  in  the  figure.  Qualitatively  the  measurements  substantiate  the 
main  features  of  the  acoustic  solution. 
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For  a quantitative  comparison  of  measured  and  predicted  overpres- 
sure signatures  the  measured  results  of  Fig.  9 have  been  reproduced  in  Fig. 

10  and  compared  directly  with  the  predicted  signatures.  Even  for  these 
extremely  strong  acoustic  waves  (amplitudes  of  about  5 kN/m?  or  one- twentieth 
of  an  atmosphere)  the  predicted  and  measured  profiles  are  in  good  agreement. 

From  experimental  data  such  as  those  shown  in  Fig.  9 it  was  found  that  the 
acoustic  analysis  was  still  valid  when  the  pressure  difference  across  the 
diaphragm  (Ap,p)  was  as  high  as  25  KN/m?,  giving  wave  amplitudes  of  approximately 
10  kN/mR.  Note  that  one  of  the  main  differences  between  measured  and  pre- 
dicted signatures  is  that  the  rapid  changes  in  overpressure  in  each  experimental 
profile  are  not  instantaneous  (see  Fig.  9)*  The  long  rise  time  of  about  0.5  ms 
associated  with  each  of  these  shocks  is  due  to  poor  diaphragm  breakage  result- 
ing from  the  relatively  small  pressure  difference  across  the  diaphragm. 


3-  PYRAMIDAL  SHOCK  TUBE 
3.1  General  Solution 


The  pyramidal  shock  tube,  which  is  depicted  in  Fig.  2b,  consists 
essentially  of  both  a pyramidal  driver  and  channel  having  the  same  divergence 
angle  and  joined  at  a common-area  station  where  a suitable  diaphragm  can  be 
installed.  This  diaphragm  separates  a normally  atmospheric  pressure  channel 
gas  from  a slightly  higher  pressure  driver  gas.  For  generality  of  the 
analysis,  which  follows  that  in  Chapter  2,  let  the  initially  quiescent  driver 
and  channel  gases  be  different  and  also  have  different  temperatures.  Before 
the  diaphragm  is  broken  the  appropriate  initial  conditions  can  be  summarized 
mathematically  as  follows. 


Channel  ( r > r ) : 

o 

n 

rH 

CL 

(3.1) 

II 

o 

(3.2) 

Driver  (o  < r < r ) : 
^ o' 

Ap2  = APq 

(3.3) 

Ai^  = 0 

(3.4) 

The  respective  symbols  Ap,  Au,  and  ApQ  denote  overpressure,  particle  velocity, 
and  pressure  difference  across  the  diaphragm,  whereas  the  respective  subscri- 
pts 1 and  2 refer  the  appropriate  symbols  to  the  channel  and  driver  gases . 

The  radial  distance  r is  measured  along  the  shock  tube  starting  from  the  apex 
of  the  pyramid  as  shown  in  Fig.  2b,  and  the  diaphragm  is  located  at  a distance 
rQ,  thereby  making  the  driver  length  equal  to  rQ. 

For  the  preceding  initial  conditions  the  ensuing  wave  motion  in  the 
pyramidal  shock  tube  after  the  diaphragm  is  broken  is  rather  complex.  This 
wave  motion  can  be  depicted  conveniently  on  a time-distance  diagram,  as  shown 
in  Fig.  6,  where  the  locus  of  each  wave  front  has  been  drawn.  On  breaking  the 
diaphragm  in  the  shock  tube  the  resulting  rapidly  expanding  driver  gas  pro- 
duces a weak  shock  wave  (g.)  in  the  channel  gas.  Simultaneously  a weak  rare- 
faction wave  (f^)  moves  into  the  driver  gas  and  eventually  reaches  and  re- 
flects from  the  cLriver  apex.  This  reflected  wave  (h.)  propagates  through 
the  driver  gas  and  interacts  with  the  contact  surface  or  interface  of  the  driver 
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and  channel  gases.  Because  the  contact- surface  displacement  is  normally 
negligible  in  c empoison  with  the  driver  length,  the  contact- surface  path  is 
shown  in  Fig.  6 sl^ly  as  a vertical  dashed  line.  Owing  to  the  change  in 
specific  impedance  (product  of  the  density  and  sound  speed)  across  the  contact 
surface  the  reflected  wave  (h_)  is  partially  transmitted  to  the  channel 
(gg)  and  partially  reflected  hack  into  the  driver  (f0).  The  subsequent 
oscillatory  wave  motion  between  the  contact  surface  and  driver  apex  creates 
a sequence  of  waves  in  the  channel.  The  integrated  result  of  all  of  these 
individual  waves  (g^,  i = 1,  2,...,oo)  yields  the  total  wave  in  the  channel. 

For  a low-pressure-ratio  pyramidal  shock  tube  for  which  the  pressure 
difference  across  the  diaphragm  Ap  is  much  less  than  the  absolute  pressure  of 
the  channel  gas,  the  waves  in  the  driver  and  channel  are  relatively  weak.  The 
wave  motion  can  consequently  be  described  adequately  by  using  well-known 
acoustic  theory.  The  one-dimensional  spherical  wave  equations  which  govern 
the  acoustic  wave  motion  in  the  channel  and  driver  gases  are  given  below. 

Channel  ( r > r ) : 
o 


Driver  (0  < r < r ) : 
o 


^(r^) 


2 

*1 


^(r^) 

"77“ 


(3.5) 


d2(H>2) 

""77" 
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a2 


<32(r<t>2) 

"“57" 


(3.6) 


The  respective  symbols  $,  a and  t denote  the  total  velocity  potential,  sound 
speed  and  time. 

For  the  analysio  it  is  convenient  to  express  the  total  velocity  po- 
tential of  the  channel  gas  (<0  in  its  more  basic  form  consisting  of  the  sum 
of  velocity  potentials  of  individual  waves  in  the  channel.  This  is  simil- 
arly so  for  the  total  velocity  potential  of  the  driver  gas  (<t>0).  The  new 
forms,  and  the  related  expressions  for  overpressure  and  particle  velocity, 
are  all  summarized  below. 


Channel  (r  > rQ): 
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Api  = -pi  3T 
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Driver  (0  < r < r ) : 
o 

Ap  t t V1  f 

»2  = - + ± ^ [q(ri)  H{t|-2i+l}  + hj(p)  H{p-2i+l} 

2 i-1 

<3<t> 

Ap2  = - P2 

d#,., 

Au2  = zr 


(3.10) 

(3.11) 

(3.12) 


In  the  above  expressions:  p denotes  density;  the  respective  nondimensional 
variables  i,  T|  and  6 equal  (a,t  - r + r )&  /a~r  , (a?t+  r)/r  and  (a  t - r)/r  ; 
g’(l)  , f ^ ( rj)  and  h!(p)  denote'the  velocity  potentials  of  the0^*1  disturbance? 
(see  fig'i  6)  and  tfee  prime  (')  denotes  differentiation  of  the  variable  with 
respect  to  the  argument  given  in  the  following  brackets;  and  H{T^]  denotes  the 
unit  step  function  which  equals  zero  prior  to  the  arrival  of  the  it*1  wave 
(Tj<  0)  and  equals  unity  after  the  arrival  of  its  wave  front  (Ti>  0)  . The 
form  of  the  total  velocity  potentials  4>-j_  and  4>2  (Eqs.  3.7  and  3-10)  are 
quite  general  but  not  arbitrary,  as  they  have  been  chosen  to  satisfy  both 
their  respective  wave  equations  (Eqs.  3-5  and  3.6)  and  their  associated  ini- 
tial conditions  (Eqs.  3*1  and  3*2  and  Eqs.  3*3  and  3 • As  an  example  of 
satisfying  the  initial  conditions  the  term  -Ap0  t/p2  in  Eq.  3-10  accounts  for 
the  pressure  difference  across  the  diaphragm  (Ap0)  for  time  t less  than  zero. 


The  presently  known  velocity  potentials  g!  (l),  f J ( tj)  and  h!  ((3)  can 
now  be  determined  by  evaluating  the  effects  of  two  different  boundaries  on  the 
wave  motion.  The  first  boundary  to  be  considered  is  the  driver  apex.  At 
this  stationary  boundary  the  particle  velocity  of  the  driver  gas  must  be 
zero  for  all  time.  By  setting  the  particle  velocity  of  the  driver  gas  Au^, 
(Eqs.  3.10  and  3*12)  equal  to  zero  at  this  boundary  where  distance  r is 
allowed  to  approach  zero  the  following  intermediate  result  can  be  obtained. 


lim 
r->  0 


£a 


YYfJU)  Hh  -21+1}  - h£(p)  H{p-2i+l) 


i=l 

00 


^2  ■2i+1^  +hi(p)  H(3-2i+i] 


i-1 


= 0 


(3.13) 


In  the  limit  when  r equals  zero  the  following  result  can  be  obtained. 


^ q(t|)  H(t)  -2i+l)  + h’(P)  H{p-2i+l) 

i-1 


= 0 


(3.1*) 


Owing  to  the  mathematical  structure  of  Eq.  3*1^  this  expression  will  in 
general  be  identically  zero  for  all  time  if  and  only  if  certain  terms  having 
equivalent  step  functions  cancel  exactly.  As  ii  and  P both  equal  at 
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the  boundary  r equal  to  zero,  the  lengthy  Eq.  3.1*+  can  consequently  be 
rewritten  in  the  following  equivalent  and  more  convenient  form. 


h!(p)  = -f'(Tj)  i=l,  2,...,oo  (3.15) 

This  result  illustrates  that  the  particle  velocity  of  the  ith  wave  having 
a velocity  potential  f * (tj)  is  countered  exactly  at  the  driver  apex  by  the 
particle  velocity  of  the  corresponding  reflected  wave  having  a velocity 
potential  h^(P),  thereby  leaving  no  net  particle  velocity  at  the  stationary 
boundary. 


The  second  boundary  to  be  considered  is  the  interface  of  the 
driver  and  channel  gases.  At  the  contact  surface  the  overpressure  and  par- 
ticle velocity  are  taken  to  be  continuous  after  the  diaphragm  is  removed, 
and  the  effects  of  diffusion,  heat  transfer  and  turbulent  mixing  are  neg- 
lected. The  matching  of  the  driver  and  channel  gas  overpressures  and 
particle  velocities  at  the  contact  surface  can  be  difficult  mathematically 
because  the  contact  surface  is  not  stationary  nor  is  its  motion  known  apriori. 
To  circumvent  this  difficulty  the  effects  of  the  contact- surface  motion  axe 
assumed  to  be  negligible  and  the  matching  procedure  is  done  at  the  dia- 
phragm station  - a fixed  or  stationary  location  r . 


For  the  first  step  in  the  matching  procedure  at  the  diaphragm 
station  (r=rQ)  the  overpressure  of  the  channel  gas  Ap,  (Eqs.  3-7  and  3.8) 
and  the  driver  gas  Ap2  (Eqs.  3*10  and  3 .11)  are  equated,  thereby  yield- 
ing the  following  intermediate  result. 


00 

iZk 

i=l 


(5)  H[|-2H2) 


■ r< A> 


+ P, 


^(l)  Hh-2i+l)  + h£0)  H[P-2i+l) 
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(3.16) 


For  the  boundary  in  question  for  which  r equals  rQ,  5 equals  a^t/r  , t) 
equals  (a^t  + r )/r  and  P equals  (a^t  - r )/r  . Certain  terms  in°Eq. 

3.l6  consequentJy  h&ve  equivalent  step  functions,  which  can  be  grouped 
accordingly,  and  expressed  in  the  following  equivalent  and  more  convenient 
form. 


P1g^(!)=  -Apo  ro/a2  + pg  f”(r,)  i-1 


(3.17) 


Px  gj(e)  = P2  hj^O)  + P2  f"(n) 


i=2,  3,...,« 


(3.18) 


The  first  expression  (Eq.  3*17)  shows  that  at  the  contact  surface  the  over- 
pressure of  the  first  (shock)  wave  in  the  chaumel  (g^)  matches  exactly  the 
sum  of  the  initial  overpressure  in  the  driver  (£pQ)  and  the  overpressure  of 
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the  first  (rarefaction)  wave  in  the  driver  (f-^) . The  second  expression 
(Eq.  3 .18)  shows  that,  for  the  i^h  wave  interaction  at  the  contact 
surface  (see  Fig.  6),  the  sum  of  the  overpressures  of  the  incident  wave 
(lu_^)  and  reflected  wave  (f^)  matches  exactly  the  overpressure  of  the 
transmitted  wave  (g^). 

Equation  3*17  and  3*18  can  be  integrated  with  respect  to  time  to 
yield  the  following  useful  results. 


P1  gl^  = " AP0  ro  t + p2  1=1  (3.19) 

Px  «i(0  = P2  hi_iW  + P2  i=2,  3,...,»  (3.20) 


The  constants  of  integration  are  each  equal  to  zero  because  the  velocity 
potentials  g^(0>  fi(ij)  and  h?  (£)  are  each  equal  to  zero  at  their  respective 
wave  fronts.  Note  that  Eqs.  13.19  and  3-20  illustrate  that  the  product  of 
density  and  velocity  potential  (p<i>)  is  continuous  at  the  contact  surface 
after  the  diaphragm  has  been  removed,  just  like  the  overpressure  (Ap)  and 
particle  velocity  (Au) . 


The  second  and  final  step  of  the  matching  procedure  at  the  dia- 
phragm station  (r=r  ) is  to  equate  the  particle  velocity  of  the  channel  gas 
Au^  (Eqs.  3.7  and  3°9)  with  that  of  the  driver  gas  Au£  (Eqs.  3.10  and  3.12), 
thereby  giving  the  following  intermediate  result . 
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00 

^ fhj(p)  + h'(P))  H{P-2i+l]  - (f£(ij)  - fj(ro)  H(T,-2i+l] 
i=l 


(3.21) 


Like  the  previous  result  for  the  overpressure  (Eq.  3»l6)  certain  terms  in 
Eq.  3*21  have  equivalent  step  functions,  which  can  be  grouped  accordingly, 
and  expressed  in  the  following  equivalent  and  more  useful  form. 


^ gJ(S)  + 6{(l)  = -f£(l)  + f’(Tl) 


i=l 


^ 6^(0  + gj(l)  = h^^P)  + h'^P)  - + q(Tl) 


i' 

i=2,  3, 


(3.22) 

(3.23) 


The  first  expression  (Eq.  3.22)  shows  that  at  the  contact  surface  the  particle 
velocity  of  the  first  (shock)  wave  in  the  channel  (g-^)  matches  that  of  the 
first  (rarefaction)  wave  in  the  driver  (f^).  The  second  expression  Eq.  3.23 
shows  that,  for  the  i™1  wave  interaction  at  the  contact  surface  (see  Fig.  6), 


19 


the  sum  of  the  particle  velocities  of  the  incident  wave  (hi.p)  and  re- 
flected wave  (fi)  matches  exactly  the  particle  velocity  of  the  trans- 
mitted wave  (g^). 

There  is  actually  a third  boundary  which  is  the  large  end  of 
the  pyramidal  channel.  However,  this  boundary  need  not  be  considered  here 
as  it  has  already  been  illicitly  assumed  that  the  channel  is  either  in- 
finitely long  or  terminated  by  a perfect  reflection  eliminator  such  that 
no  reflected  or  other  waves  arise  from  this  boundary. 


Now  that  the  matching  procedure  has  been  completed  the  resulting 
differential  equations  for  gi(|),  f (ij)  and  hi(p)  can  be  cast  into  a more 
amendable  form.  From  Eqs.  3.17,  3.19  and  3.22  one  can  obtain  differential 
equations  for  gpd)  and  f i(tj)  . In  a similar  manner,  from  Eqs.  3.l8.  3.20 
and  3*23,  differential  equations  can  be  obtained  for  gj.(?)  and  Mr]),  where 
i is  greater  than  unity.  These  results  associated  with  the  boundary  r 
equal  to  rQ  and  the  p- evious  result  (Eq.  3.15)  for  the  other  boundary  r 
equal  to  zero  are  summarized  below. 


g£(D 


fj(i) 


a (p  -p  ) 

+ — - _ o'  ( f ) 

alPl+a2P2 


a-,(p2"pi) 

+ ■■■  -I--  f*; (tj) 

Vl+a2p2  1 


( Ap  rc  a, 

— - — — ( |-1) 

a1p1+a2p2  a2 


2a2p2 


£ c 1 >.'1  /g\ 

a1P1+a2p2  a2  i-lVP; 


alpl+a2P2  *2 
aiPl'a2P2  „ , . 

Vl+a2P2  i-l( 


i=l 


i=2,  3,  • • . ,c 


i=l 


(3.24) 


1l(P2-pL) 

- h! 


(3.25) 


alpl+a2p2 


i-l<« 


i=2,  3,. 


C» 


h|(P)  - -f’(Tl) 


i-1,  2,, 


(3-26) 


The  above  set  of  differential  equations  for  g^l),  fi(*))  and  hj^P)  can 
be  uncoupled  and  subsequently  solved  to  give  explicit  expressions  for 
6i(0  , fih)  and  h^p)  . 


A convenient  initial  step  in  solving  the  differential  equations 
is  to  find  a recurrence  relationship  for  f'(t|)  . Using  Eq.  3.26  we  can 
rewrite  Eq.  3-25  in  the  following  form. 
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The  nondimensional  parameters  b and  l are  equal  to  agPs/aqPi  and 
(aiPi+a2P2)/al(P2“Pl)  respectively.  It  should  be  noted  that  f^_i(T)-2) 
and  fj,(rj-2)  in  Eq.  3-27  are  expressed  correctly,  as  the  variable  tj-2 
in  the  brackets  (instead  of  simply  r\)  is  required  to  ensure  that  the  func- 
tion fji  (t)-2)  on  the  right-hand  side  of  the  equation  is  properly  phased  in 
time  with  the  function  f^ ( tj)  on  the  left-hand  side.  After  multiplying  the 
ith  differential  equation  for  fi(l])  (Eq.  3*27)  by  exp  [(i)-2i+l )/£]  and 
partially  integrating  with  respect  to  time  the  following  recurrence  relation 
for  f{(n)  can  be  obtained. 


f;(n)  = 


Ap  r l 
*o  o 

a2P2(B  + 1) 


(i1-l)/i 


exp 


(-  ^ J [(l+B+i)  y + h]  exp(y)  dy  i=l 


J o 

* -k  ^ ) J 


(H-2i+l)/i  (3,28) 
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From  these  results  the  i velocity  potential  can  be  determined,  and  it  is 
given  below. 


fi<">  - i 1)  X 0-i)  ( *r  (3-29) 

J-l 

The  notation  (n)  denotes  the  well-known  binomial  coefficients,  which  can  be 
expressed  as  m'./n'.(m-n) '. , where  the  symbol  '.  denotes  the  factorial  function. 
For  convenience  and  illustrative  purposes  the  first  five  velocity  potentials 
f/(ij)  are  given  in  expanded  form  in  Appendix  A,  along  with  the  corresponding 
first  four  expanded  functions  I'(ij-2i+l).  The  functions  I.(i)-2i+l)  can  be 
expressed  in  the  following  recurrence  form.  J 

r u/t, 

2xp(-u/i)  J [(1+B+iV  +b]  exp(y)  dy  j=l 

\/l  <3-3°) 

J=2,  3 °° 
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p “V* 

exp(-u/£)J  exp(y)  dy 
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Alternatively  the  j function  I'(ij-2i+l)  cm  be  expressed  in  the  following 
more  elegant  form.  d 


/u/jt  p u>/ 1 

...j...  J [(1+b+j0)  y + a]  exp(y)  dy  (3-31) 

In  either  case  the  final  results  for  the  function  I'(tj-21+1)  can  be 
determined,  and  it  is  given  below.  J 

/ \ A k(l+a+/)-B  r \ J-k 

lj(u>)  * (1-Hs+i)  ( j - ij  + ■ + ^ '( j.k)  ■ — ( J -J  exp(-w/ /)  (3-32) 
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The  solution  for  the  wave  motion  in  the  pyramidal  shook  txibe  utilizing 
different  driver  and  channel  gases  has  now  been  obtained,  although  further  work 
is  required  to  get  explicit  expressions  for  the  overpressure,  particle  velocity 
and  particle  path  of  the  driver  and  channel  gases.  The  result  for  the  i*'*1 
velocity  potential  f£(ti)  as  given  by  Eqs.  3*29  and  5»32  can  now  be  used  to 
obtain  the  i"th  velocity  potentials  for  h£(P)  and  g£(|).  The  appropriate 
expressions  for  h^(0)  and  g^(|)  which  stem  from  Eqs.  3.15,  3.19  and  3-20,  are 
given  below. 
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These  results  for  g£(|),  fi(T))  and  hj[(B)  can  be  substituted  into  the  expressions 
for  the  total  velocity  potentials  0^  and  0o  (Eqs.  3.7  and  3»10).  Then  the  two 
respective  expressions  for  the  overpressure  and  particle  velocity  of  the  channel 
gas  follow  from  Eqs.  3.8  and  3»9i  and  those  for  the  driver  gas  follow  from 
Eqs.  3.11  and  3.12.  Although  the  solution  is  rather  complex  and  lengthy,  the 
overpressure  and  particle  velocity  signatures  for  both  the  channel  and  driver 
gases  can  be  obtained  fairly  easily  by  using  a digital  computer. 


It  can  be  seen  from  the  solution  of  the  wave  motion  in  the  pyramidal 
shock  tube  that  both  the  ratios  of  the  driver  and  channel  gas  sound  speeds 
( ®2/ ai)  and  densities  (po/p^)  play  important  roles  not  only  separately  but 
also  in  a combined  form  (e.g.,  see  Eqs.  3.29>  3»32  and  3«3M«  The  important 
confcination  agp^a^p^  has  been  denoted  by  the  symbol  et.  This  nondimen sional 
parameter  car  be  expressed  in  the  following  alternate  forms  by  using  the 
equation  of  state  (p  - pRT)  and  the  isentropic  sound-speed  relation  (a^  = pp/p  ~ 
7Rot/M). 
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The  respective  symbols  y,  p,  R,  T,  Rc  and  M denote  the  specific  heat  ratio, 
pressure,  gas  constant,  temperature , universal  gas  constant  (R0  = MR)  and 
molecular  weight.  It  is  worth  noting  that  the  specific  impedance  ratio  b 
equals  unity  not  only  for  the  simplest  case  of  identical  driver  and  channel 
gases  having  the  same  tenperatures  but  also  for  the  specific  case  of  different 
driver  and  channel  gases  provided  that  M^/Ti  equals  72^/^2*  Furthermore, 
the  specific  impedance  ratio  ■ equals  (72^/71*tL/  ' ^ Tor  the  very  important 
case  of  different  driver  and  channel  gases  having  equal  temperatures  and  b 
equals  (Tq/T^)V2  for  the  not  so  practicable  case  of  identical  driver  and 
channel  gases  having  different  temperatures. 

A set  of  values  of  the  specific  impedance  ratio  b was  compiled  for 
certain  different  combinations  of  driver  and  channel  gases  having  identical 
tenperatures,  and  they  are  given  in  Table  2.  The  smallest  values  of  b occur 
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when  a light  driver  gas  such  as  hydrogen  (H2)  or  helium  (He)  is  used  in 
conjunction  with  a heavy  channel  gas  such  as  dichlorodifluoromethane  (CCI2F2) , 
sulfur-hexafluoride  (SFg)  or  octofluorocyclobutane  (Cl+Fg).  On  the  other  hand 
the  largest  values  of  b occur  when  the  driver  gas  is  heavy  or  has  a large 
molecular  weight  and  the  channel  gas  has  a low  molecular  weight.  For  nearly 
equivalent  molecular  weights  of  the  driver  and  channel  gases  the  specific 
inpedance  ratio  is  always  nearly  equal  to  unity. 


Another  inport ant  combination  of  sound  speed  and  density  r atios, 
(ajP!  + a2P2)/ai(P2  - Pi)  or  (l  + a2P2/aiPi)/(P2/Pl  - l)>  has  been  denoted 
by  the  symbol  i.  This  nondimensional  parameter  can  be  expressed  in  the 
following  alternate  forms. 


1 + 


l = 


alpl 


' (L  , ' L 

1 B 

r\  v 


1 + B _ 
1 2 . 


*>T1 

“lT2 


112 


(3-36) 


- 1 


The  parameter  l,  which  plays  the  role  essentially  of  a nondimensional  time 
constant  for  the  wave  motion  (e.g.,  see  Eq.  3*32),  can  take  on  values  that  range 
from  positive  to  negative  infinity.  A set  of  values  of  the  parameter  l was 
compiled  for  certain  interesting  combinations  of  driver  and  channel  gases 
having  identical  temperatures,  and  these  results  are  shown  for  interest  in 
Table  3. 


The  wave  that  propagates  in  the  channel  gas  is  one  of  the  most 
important  features  of  the  wave  motion  in  the  pyramidal  shock  tube.  For  this 
reason  the  overpressure  and  particle  velocity  signature  of  this  wave  will  be 
explored.  The  mathematical  expression  for  the  overpressure  signature,  resulting 
from  Eqs.  3.7,  3.8}  3.29,  3-32,  3.34,  3.35  and  3-36,  can  be  expressed  as  given 
below. 
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For  convenience  and  illustration  the  first  five  functions  of  both  gi( |)  and 
Jj( |-2i+2)  are  given  expanded  form  in  Appendix  B.  Note  that  the  function  g^(|) 


23 


i 


I 


I 

1 

i 


used  in  Eqs.  3-37  and  3-38  has  been  redefined  slightly,  differing  from  the 
previous  one  by  a factor  of  /^>0r02/a2Pi(z+l) . 

Fifteen  nondimensional  overpressure  signatures  for  the  channel  of  a 
pyramidal  shock  tube,  as  computed  from  Eqs.  3.37,  3.38  and  3.39*  are  displayed 
in  Fig.  11  for  fifteen  different  combinations  of  driver  and  channel  gases  having 
the  same  temperature.  For  each  signature  the  driver  and  channel  gases  are 
depicted  in  the  insert.  The  fifteen  signatures  have  been  arranged  in  a sequential 
order  based  on  the  specific  impedance  ratio  b (shown  in  the  figure).  Consequently, 
the  combination  of  a light  driver  gas  and  a heavy  channel  gas  (z  < 1)  appear  first 
in  the  sequence  (Fig.  11a  to  lie),  and  the  combination  of  a heavy  driver  gas  and 
a light  channel  gas  (z  > l)  appears  last  (Fig.  llg  to  llo).  It  can  be  seen  that 
an  overpressure  signature  corresponding  to  a a value  that  is  small  (z  « l)  or 
large  (a  » l)  differs  markedly  from  an  ideal  N-shaped  profile  for  a equal  to 
unity  (Fig.  Ilf).  However,  as  both  an  initially  small  or  large  value  of  b 
approaches  and  equals  unity  the  differences  decrease  and  disappear.  For 
equivalent-temperature  driver  and  channel  gases,  therefore,  only  identical  driver 
and  channel  gases  can  produce  an  N-shaped  overpressure  signature  in  the  channel 
of  a pyramidal  shock  tube.  This  conclusion  can  be  deduced  directly  from  Eqs. 

3*37,  3.38  and  3.39,  although  it  is  not  obvious. 

Five  nondimensional  overpressure  signatures  for  the  channel  of  a 
pyramidal  shock  tube,  as  computed  from  Eqs.  3*37,  3.38  and  3.39,  are  shown  in 
Fig.  12  for  five  different  cases  of  identical  driver  and  channel  gases  having 
different  temperatures.  The  five  driver  gas  to  channel  gas  temper ature  ratios 
T2/T1  are  9,  2,  1,  l/2  and  l/9,  as  shown  in  the  figure,  corresponding  to  a 
values  of  0.33,  0.71,  1.00,  1.41  and  3.00,  respectively.  The  signatures  in  Fig. 
12a  and  12b  are  similar  to  those  for  a combination  of  a light  driver  gas  and  a 
different  equivalent-temperature  channel  gas  (see  Fig.  lid  and  lie),  because 
heating  the  driver  gas  or  alternatively  cooling  the  channel  gas  makes  the  driver 
gas  light  in  comparison  to  the  channel  gas.  On  the  other  hand,  because  cooling 
the  driver  gas  makes  it  heavier  than  the  channel  gas  or  alternatively  heating 
the  channel  gas  makes  it  light  than  the  driver  gas,  the  signatures  in  Fig.  12d 
and  12e  resemble  those  for  a combination  of  a heavy  driver  gas  and  a different 
equl valent-temperature  channel  gas  (see  Fig.  Hi  and  11m).  For  the  case  of 
identical  driver  and  channel  gases  having  the  same  temperature  the  signature 
is,  of  course,  N shaped  as  shown  in  Fig.  12c. 

It  is  of  interest  to  consider  the  overpressure  signature  in  the  channel 
of  a pyramidal  shock  tube  for  the  very  special  case  when  different  driver  and 
channel  gases  are  used  but  they  have  equivalent  densities.  For  this  special  case 
the  parameter  / (Eq.  3*36)  is  infinite,  and  the  parameter  s (Eq.  3«35).  equals 
&2./&1  which  is  not  arbitrary  but  can  be  shown  to  be  equal  to  (72/71)^' Further- 
more, to  achieve  such  a special  case  of  identical  densities  for  different  driver 
and  channel  gases  the  temperature  ratio  Tg/T^  must  be  equal  to  Five  over- 

pressure signatures  corresponding  to  values  of  ■,  ag/a^  or  (72/7])1/2  of  0.89, 

0.8l,  1.00,  1.13  and  1.23,  as  computed  from  Eqs.  3.37,  3.38  and  3*39,  are 
displayed  in  Fig.  13.  Note  that  specific  sets  of  values  for  72  and  7^  could 
be  1.67  and  1.10,  1.40  yid  1.10,  1.10  and  1.40,  and  1.10  and  I.67  to  give  values 
of  ■ , ag/a^  or  (72/71)1/2  equal  to  9*89,  0.8l,  1.13  and  1.23,  respectively.  The 
most  interesting  feature  of  these  signatures  (Fig.  13)  is  that  the  variation  in 
overpressure  is  linear  between  successive  sudden  changes  in  overpressure  (shocks). 
Such  a linear  variation  occurs  only  when  the  parameter  /is  infinite  or  when 
the  driver  and  channel  densities  are  equal.  For  this  special  case  of  equivalent- 
density  driver  and  channel  gases  a single  N-wave  results  only  when  the  sound 
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speeds  are  also  equal  (see  Fig.  13c).  Hence,  we  note  once  again  that  a single 
N-shaped  overpressure  signature  is  possible  with  a pyramidal  shock  tube  only 
if  the  driver  and  channel  are  identical  and  have  the  same  temperature. 

The  overpressure  signature  predicted  by  Eqs.  3*37,  3*38  and  3. 39, 
some  of  which  have  been  shown  in  Figs.  11,  12  and  13,  are  invariant  with 
radial  distance  (wave  form  does  not  change  shape)  but  attenuates  inversely 
with  distance  (acoustic  l/r  decay  law).  The  amplitude  of  the  first  shock 
is  given  by  Ap0r0/(B+l)r . Consequently,  for  a given  pressure  difference 
across  the  diaphragm  (A?0),  a combination  of  a light  driver  gas  and  a heavy 
channel  gas  having  a small  specific  impedance  ratio  (b)  produces  a relatively 
strong  wave  in  the  channel.  However,  the  duration  of  this  wave  is  generally 
very  short  because  the  time  interval  between  successive  shocks,  is 

small  owing  to  the  high  sound  speed  ( ) of  a light  driver  gas.  On  the 
other  hand,  the  relatively  low- amplitude  wave  for  a combination  of  a heavy 
driver  gas  and  a light  channel  gas  can  have  a very  long  duration. 

Although  the  overpressure  signature  described  by  Eqs.  3*37,  3-38 
and  3.39  is  tedious  to  calculate  in  detail,  a simple  expression  can  be  found 
to  give  the  sudden  change  in  overpressure  associated  with  each  shock  in  the 
signature.  Firstly,  since  | equals  2i-2  at  the  wave  front  of  the  ith  distur- 
bance whgre  the  shock  is  located,  find  fromtiEq.  3*39  that  Jo(0)  is  equal  to 
unity,  J^(0)  eguals  h,  and  a^l  subsequent  jj(0)  are  zero.  Secondly,  from 
Eq.  3*38  for  gi(i),  obtain  gi(0)  equal  to  unity  and  subsequent  gV(2i-2)  equal 
to  2B(B-l)i”2/tB+l)i“l.  Finally,  find  from  Eq.  3-37  the  desired  result  for 
the  amplitude  of  the  i^*1  shock,  which  is  given  below. 
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Relative  to  the  first  shock  having  an  amplitude  of  <£p0r0/(a+l)r,  the  amplitude 
of  the  second  one  is  2b/(b+1).  Therefore,  depending  on  the  value  of  the 
specific  impedance  ratio  a,  the  second  shock  can  be  smaller  (a  < l),  equal 
(b  = l),  or  larger  (b  > l)  in  anplitude  than  the  first  one.  However,  regardless 
of  the  value  of  b,  the  amplitude  (absolute)  of  the  third  and  subsequent  shocks 
always  diminish,  being  smaller  than  its  predecessor  by  the  factor  (b-1)/(b+1). 
(see  the  overpressure  profiles  given  in  Figs.  11,  12  and  13.) 

The  mathematical  expression  for  the  particle-velocity  signature  of 
the  wave  propagating  in  the  channel  gas  can  be  determined  directly  from  Eqs. 

3*7,  3*9,  3*29»  3.32,  3 *3^ > 3.35  and  3*36.  Alternatively  and  more  simply  it 
can  be  derived  by  using  the  following  acoustic  expression  which  relates  the 
particle  velocity  to  the  known  overpressure  (Eqs.  3*37,  3-38  and  3*39)  of  a 
spherical  wave  moving  away  from  the  origin. 

^ ^ + pr  f ^(y-r/a)  *y  (3.M) 

-CO 


Ap  r 
~2_o_o 

(■+!)* 


'Vo  f 2b 

(a+l)r  V b+1 


g-1 

B+1 


i-2 


25 


In  either  case  the  final  results  for  the  particle-velocity  signature  are  given 
below. 
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Note  that  the  i function  gi(|)  and  the  associated  function  Jj(w)  are 
not  given  here  as  they  have  been  defined  previously  by  Eqs.  3.38  and  3-39, 
respectively. 


The  particle-velocity  profile  is  not  invariant  with  radial  distance 
like  the  overpressure  profile,  because  it  consists  of  two  parts  or  two  sums 
of  terms  which  have  different  dependences  on  radial  distance.  The  part  of 
the  particle-velocity  expression  containing  the  sum  of  the  g^(|)  terms  has 
a dependence  on  radius  of  l/r2.  This  part  of  the  particle  velocity  conse- 
quently diminishes  relatively  rapidly  with  distance  and  is  therefore 
significant  only  at  small  radii  or  in  the  near-field.  The  other  part  of 
the  particle-velocity  expression  has  a weaker  dependence  on  radius  of  l/r 
and  thus  becomes  dominant  at  large  radii  or  in  the  far-field.  It  can  be 
seen  by  conparing  Eqs.  3*37  and  3.42  for  the  overpressure  and  particle 
velocity  respectively  that  the  changing  particle-velocity  profile  tends  with 
increasing  distance  to  take  on  the  shape  of  the  overpressure  profile.  For 
the  far-field  one  consequently  finds  that  Aup  is  singly  equal  to  £pi/&±p±, 
a well-known  acoustic  result.  The  parameter  in  Eq.  3*42  that  is  important 
as  a rough  indicator  when  this  far-field  result  can  be  applied  instead  of 
the  more  complete  Eq.  3.42  is  a^r^agr.  From  detailed  calculations  of 
particle-velocity  signatures  it  can  be  shown  that  when  a^r^agr  is  less 
than  about  l/4,  or  when  radial  distances  r are  greater  than  4a]r0/a2,  the 
far-field  approximation  (Au^  = APi/aiPq)  is  applicable.  It  is  worth  noting 
that,  as  the  wave  length  X between  successive  shocks  in  the  overpressure  or 
particle-velocity  signature  is  singly  2r0a^/a2s  the  far-field  result  applies 
for  radial  distances  r greater  than  about  2X. 

Nondimensional  particle-velocity  signatures  at  five  different  radial 
distances  r of  rOJ  2rQ,  4rQ,  8r0  and  oo  in  the  channel  of  a pyramidal  shock 
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tube  have  been  computed  from  Eqs.  3.42,  3.43  and  3.44  for  four  different 
combinations  of  driver  and  channel  gases,  and  these  four  sets  of  results 
are  given  in  Figs.  l4,  15,  16  and  17.  The  driver  and  channel  gases,  which 
were  taken  to  have  the  same  temperature , are  shown  at  the  top  of  each 
figure.  From  the  results  in  each  figure,  one  can  see  that  the  change  in 
shape  of  etch  signature  with  increasing  distance  from  the  diaphragm  station 
(r=r0)  is  noticeable  only  in  profiles  corresponding  to  r/r0  values  of  less 
than  about  2~K  or  2r0a^/a2,  or  values  of  0.7ro,  2rQ,  4rQ  and  l8r0  for  the 
He- Air,  Air-Air  (identical  driver  and  channel  gases),  SFg-Air  and  Cl+Fg-He 
driver -to- channel  gas  combinations  respectively.  This  waveform  distortion 
persists  to  large  radii  for  the  case  of  a high  ap/a^  value  which  corresponds 
to  a combination  of  a heavy  driver  gas  and  a light  channel  gas. 

It  is  worth  mentioning  that  a simple  expression  can  be  derived  to 
give  the  sudden  change  in  particle  velocity  associated  with  each  shock  in 
the  signature.  The  derivation  is  similar  to  that  for  the  sudden  change  in 
overpressure  associated  with  each  shock,  which  yielded  Eq.  3«40.  The 
desired  result  for  the  amplitude  of  the  i^*1  shock  in  the  particle-velocity 
signature  is  given  below. 
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These  results  are  very  similar  to  those  of  Eq.  3.40. 
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(3.45) 


The  motion  of  the  channel  gas  and  contact  surface,  owing  to  the 
wave  motion  in  the  pyramidal  shock  tube,  can  be  determined  by  using  the 
previously  derived  expression  for  the  particle  velocity  (Eqs.  3.42,  3.43  and 
3.44).  As  AU]_  is  equal  to  dr/dt  exactly  or  a£<i(r/]&)/d!  to  first  order  (same 
order  of  approximation  as  the  acoustic  equations  approximate  the  gasdynamic 
equations),  the  particle  displacement  can  be  obtained  by  integration.  For 
a fluid  particle  having  an  initial  location  r*  (r*  > rQ)  the  first-order 
particle  displacement  r-r#  or  simply  Ar  can  be  shown  to  be  given  by  the 
following  results. 
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Note  that  the  i function  g£(|)  and  the  associated  j function  Jj(w)  are 
not  given  here  because  they  have  been  defined  previously  by  Eqs.  3-^3  and 
3.41*  respectively. 


The  displacement  of  a fluid  particle  as  a fraction  of  time,  including 
the  total  displacement  after  the  wave  motion  has  subsided  at  large  times  (£-><»), 
can  be  calculated  from  Eqs.  3.^6,  3.^7  and  3.^8  quite  easily  by  using  a digital 
computer.  Such  confuted  results  will  be  given  subsequently.  First,  however, 
an  inportant,  simple  expression  will  be  derived  for  the  total  displacement  f 
a fluid  particle.  Let  the  driver  gas  have  an  initial  pressure  of  p-^+ApQ, 
where  p^  is  the  initial  and  final  pressure  of  the  channel  gas  and  Ap0  is  the 
pressure  difference  across  the  diaphragm.  For  a driver  gas  having  a total 
mass  iQq  and  an  initial  total  volume  VQ,  the  initial  density  is  simply 
Also,  for  a pyramidal  driver  having  a spherically  shaped  diaphragm  of  radius 
rD  the  driver  volume  is  directly  proportional  to  rQ3.  Now  let  the  driver  gas 
expand  isentropically  (as  it  would  for  acoustic  wave  motion)  to  a final  pressure 
Pj_.  The  corresponding  final  density  would  be  nio/v,  where  the  final  volume  V 
is  proportional  to  r3.  The  following  isentropic  expression  relating  the 
initial  and  final  pressures  and  densities  can  therefore  be  applied. 
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From  this  result  the  total  displacement  of  the  contact  surface  or  interface 
of  the  driver  and  channel  gases,  which  is  r-rQ  or  singly  ArQ,  is  given  to 
first  order  by  the  following  expression. 
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This  result  for  the  contact  surface  can  now  be  generalized  for  any  fluid 
particle  in  the  channel.  For  a contact- surface  displacement  from  radius 
ro  '*;o  ro+^ro  or  sinply  Ar0,  the  volume  of  fluid  displaced  is  directly 
proportional  to  (r0+Ar0)3  - rQ3  or  r0^Ar0  to  first  order.  The  corresponding 
displacement  of  a fluid  particle  in  the  channel  from  radius  r#  to  r#+Ar  or 
sinply  Ar  produces  an  equivalent  volume  displacement  (final  pressure  and 
density  are  spatially  uniform  in  the  channel)  which  is  directly  proportional 
to  (r#+Ar)3  - r#3  Qr  r 2Ar  to  first  order.  We  therefore  obtain  the  result 
that  r02ArQ  equals  r*2Ar.  By  using  this  result  in  conjunction  with  Eq.  3*50, 
the  total  displacement  of  a fluid  particle  Ar  from  its  initial  location  r#  in 
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the  channel  is  given  by  the  following  expression. 
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The  expressions  for  the  total  displacement  of  the  contact  surface  (Eq.  3.50) 
and  a fluid  particle  in  the  channel  (Eq.  3*51)  are  important  not  only  because 
of  their  sinplicity  but  also  they  indicate  how  to  properly  normalize  the 
full  expression  for  the  particle  displacement  (see  Eq.  3*46). 


It  would  be  of  interest  to  derive  a sinple  expression  from  Eqs . 
3.46,  3-47  and  3-48  for  the  total  displacement  of  a fluid  particle  and  check 
to  see  if  this  result  is  identical  to  Eq.  3. 51.  Such  a derived  expression 
proved  elusive.  Therefore  a digital  computer  was  used  to  calculate  (Eqs. 
3.46,  3.47  and  3.48)  the  particle  displacement  at  sufficiently  large  times 
to  obtain  the  total  displacement.  In  each  case  the  computed  result  was  in 
agreement  with  the  prediction  of  Eq.  3-51.  In  addition  to  providing  a means 
of  simply  normalizing  Eq.  3-46,  therefore,  Eq.  3.51  provides  an  excellent 
check  on  the  analysis  of  the  wave  motion  in  the  pyramidal  shock  tube. 

It  is  worth  noting  here  that  the  total  displacement  of  a fluid 
particle  Ar/r0  or  ^0ro^/372Plr-#^  relatively  small,  and  it  is  also 
independent  of  both  the  specific  impedance  ratio  e and  the  parameter  i. 
Consider,  for  example,  the  total  displacement  of  the  contact  surface,  for 
which  equals  rQ.  Even  for  a relatively  large  pressure  difference  across 
the  diaphragm  (Ap0)  of  one-twentieth  of  an  atmosphere  (p^/20),  and  70  equal 
to  1.4  for  a diatomic  driver  gas,  the  total  displacement  is  only  r0/84  or 
1.2$  of  the  driver  length.  In  general  the  contact- surface  displacement  is 
small  and  a fluid-particle  displacement  at  a larger  radius  is  even  smaller. 
Consequently  the  previous  assumption  for  the  wave-motion  and  displacement 
analyses  that  the  contact- surface  displacement  was  negligible  conpared  to 
the  driver  length  was  very  reasonable. 


Ncndimensional  fluid-particle  histories  for  different  initial 
radii  r#  of  rQ,  2r0,  4rQ  and  8r0  in  the  channel  of  a pyramidal  shock  tube 
have  been  computed  for  four  different  confcinations  of  driver  and  channel 
gases,  and  these  four  sets  of  results  are  given  in  Figs.  18,  19,  20  and  21. 

The  driver  and  channel  gases,  which  were  taken  to  have  the  same  tenperatures, 
are  shown  in  each  figure.  Figures  18  and  19,  for  driver-to-channel  gas  com- 
binations of  He- Air  and  Air-Air,  respectively,  show  that  the  particle  displace- 
ment increases  and  then  decreases  to  its  final  displacement  (asymptote  indicated). 
This  behaviour  is  typical  of  the  case  of  a light  driver  gas  and  a heavy  channel 
gas  for  which  the  specific  inpedance  ratio  es  is  less  than  unity.  From  the 
results  of  Figs.  20  and  21  for  driver-to-channel  gas  combinations  of  Ci^Fg-Air 
and  C02-He  respectively,  the  particle  displacement  increases  in  a "darped" 
oscillatory  manner  to  its  final  displacement.  This  behaviour  is  typical  of 
the  case  of  a heavy  driver  gas  coupled  to  a light  channel  gas  for  which  the 
specific  impedance  ratio  b is  greater  than  unity.  In  each  of  the  Figs.  l8, 

19,  20  and  21  the  top  displacement  history  for  radius  r#  equal  to  rQ  is  that 
for  the  contact  surface.  The  contact- surface  displacement  is  the  largest, 
and  particle  displacements  at  increasing  radii  continually  decrease  and  become, 
eventually,  inconspicuous,  because  the  wave  motion  being  strongest  at  the 
diaphragm  station  becomes  weaker  with  increasing  distance  and  eventually 
vanishes. 
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The  particle-displacement  histories  for  the  previous  driver-to- 
channel  gas  combinations  of  He-Air  and  Air-Air  given  in  Figs.  18  and  19  are 
also  shown  in  Figs.  22  and  23.  In  these  latter  figures,  however,  the  dis- 
placement histories  are  superposed  on  a time- distance  diagram  which  also 
depicts  the  wave  motion.  Additionally,  displacement  histories  for  another 
driver -to- channel  gas  combination  of  SFg-Air  are  similarly  shown  for  interest 
in  Fig.  2k  . 


In  order  to  check  the  validity  of  the  acoustic  analysis  of  the  wave 
motion  .in  the  pyramidal  shock  tube,  the  UTIAS  Travelling- Wave  Horn  was 
used  to  obtain  experimental  data  which  could  be  compared  with  predicted 
results.  The  pyramidal  shock  tube  had  a total  divergence  angle  of  7.2 
degrees,  an  overall  length  of  25m  and  a driver  that  was  1.58m  long  (Refs. 

3 and  5) . The  driver  gases  used  for  the  experiments  were  helium  (He),  air, 
argon  (A),  carbon  dioxide  (CC^),  sulfur  dioxide  (SO2),  dichlorodifluoromethane 
(CCI2F2)  and  octafluorocyclobutane  (C^Fg),  which  covered  a large  range  of 
both  molecular  weights  and  sound  speeds  (see  Table  2).  For  all  experiments 
the  channel  gas  was  air,  and  the  driver  and  channel  gases  had  the  same 
initial  temperature . Overpressure  measurements  made  in  the  channel  of  the 
shock  tube  for  the  different  driver  gases  are  shown  in  the  top  part  of 
Figs.  25  to  31-  Note  that  the  driver  and  channel  gases,  pressure  difference 
across  the  diaphragm,  driver  length  and  measurement  location  are  given  above 
each  oscillogram.  Each  experimental  overpressure  signature  has  been  re- 
produced below  its  oscillogram,  where  it  is  conpared  directly  to  the  predicted 
signature.  For  each  one  of  the  seven  comparisons  the  predicted  and  measured 
profiles  are  in  good  agreement,  substantiating  the  analysis  for  the  wave 
motion  in  the  pyramidal  shock  tube.  From  experimental  data  such  as  those 
shown  in  Figs.  25  to  31  it  was  found  that  the  acoustic  analysis  was  still 
valid  when  the  pressure  difference  across  the  diaphragm  (£pQ)  was  as  high 
as  15  kN/m2. 


3 .2  Special  Solution 

It  is  of  academic  interest  to  investigate  the  solution  of  the  wave 
motion  in  a pyramidal  shock  tube  for  the  special  case  of  different  driver 
and  channel  gases  which  have  the  same  density.  For  this  special  case  of 
equivalent  densities  the  ratio  of  the  driver  and  channel  gas  sound  speeds 

and  the  specific  impedance  ratio  b (a^^/a^p^) , are  both  fixed,  being 
equal  to  the  square  root  of  the  specific  heats  ratio  ( 72/71) - Furthermore, 
the  temperature  of  the  driver  and  channel  gases  cannot  be  arbitrary,  as  the 
temperature  ratio  T2/T;l  must  be  equal  to  the  ratio  of  the  molecular  weights 
Kq/M i- 


For  this  special  case  of  equivalent  driver  and  channel  gas  densities 
the  solution  of  the  wave  motion  takes  a relatively  single  form.  Certain 
terms  vanish  from  the  differential  equations  for  the  velocity  potentials 
gi(?),  f i (tj ) and  h^(P)  (see  Eqs.  3-2^,  3-25  sued  3-26).  The  resulting  simpler 
differential  equations  are  given  below. 
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From  these  results  the  following  recurrence  relationships  for  gV(0>  f V ( tj) 
and  h^(p)  can  be  obtained. 
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These  expressions  can  be  integrated  to  obtain  recurrence  relationships  for 
gJ(|),  f^(t))  and  hJ(P),  and  also  for  gj^(i),  fj.(rj)  and  h^P) . The  results  for 
gj(|),  f^(tj)  and  h£(P)  can  be  substituted  into  the  expressions  for  the  total 
velocity  potentials  0]_  and  02  (Eqs.  3.7  and  3 .10).  Then  the  two  respective 
expressions  for  the  overpressure  and  particle  velocity  of  the  channel  gas 
follow  from  Eqs.  3.8  and  3-9>  whereas  those  for  the  driver  gas  follow  from 
Eqs.  3.11  and  3.12. 


The  mathematical  expression  for  the  overpressure  signature  of  the 
wave  moving  in  the  channel  of  a pyramidal  shock  tube  is  given  below. 
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From  this  expression  it  can  be  shown  that  the  overpressure  signature  consists 
of  a sequence  of  diminishing  anplitude  N-shaped  segments,  some  of  which  may 
be  inverted,  but  each  having  the  same  duration  of  2t0/bq.  Five  overpressure 
signatures  corresponding  to  values  of  b,  or  ( 72/71)  V2  of  O.89,  0.8l, 
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1 00,  1.13  and  1.23  are  shown  in  Fig.  13,  illustrating  the  preceding  remarks. 
Although  these  signatures  could  have  been  generated  by  using  Eq.  3.58,  they 
were  originally  confuted  from  equivalent  expressions  of  Eqs.  3.37,  3.38  and 
3.39  of  Section  3.1.  Note  that  the  signatures  have  been  discussed  previously 
in  Section  3.1. 

The  linear  variation  in  overpressure  of  the  i^*1  N- shaped  segment  of 
an  overpressure  signature  can  be  expressed  quite  singly  as  shown  below. 
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This  result,  which  can  be  derived  from  Eq.  3-58,  is  much  easier  to  use  than 
Eq.  3.58  when  constructing  an  overpressure  signature.  From  Eq.  3.59  (or 
Eq.  3.58)  one  can  show  that  the  slope  of  the  i^h  linearly  varying  over- 
pressure segment,  or  the  first  derivative  of  the  overpressure  with  respect 
to  nondimensional  time  §,  is  given  as  follows. 
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Note  that  the  absolute  magnitude  of  the  slope  diminishes  in  the  form  of  a 
geometric  progression  having  the  common  ratio  (b-1)/(b+1). 

The  particle-velocity  signature  of  the  wave  propagating  in  the 
channel  gas  of  a pyramidal  shock  tube  can  be  determined  directly  from  Eqs. 

3.9  and  3*57.  However,  it  can  also  be  obtained  from  the  overpressure 
signature  (Eq.  3-58  or  3.59)  by  using  Eq.  3.^1.  The  easiest  procedure  is  to 
use  Eqs.  3.^1  and  3.59-  Then  the  particle  velocity  associated  with  the  i^h 
segment  of  the  signature  is  given  by  the  following  expression. 
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Each  segment  of  the  signature  consists  of  a near-field  and  a far-field  term. 
The  particle-velocity  signature  therefore  tends  with  increasing  radial 
distance  to  take  on  the  same  shape  as  the  overpressure.  Generally  the  near- 
field term  is  unimportant  when  the  parameter  a^r  /agr  is  less  than  l/U,  or 
when  the  radial  distance  r is  greater  than  i*r0a2/a^  or  2A, where  A is  the 
wave  length  of  each  segment. 

The  motion  of  the  channel  gas  and  contact  surface,  owing  to  the 
wave  motion  in  the  pyramidal  shock  tube,  can  be  determined  from  Eq.  3.61  for 
the  particle  velocity.  As  Aui  is  equal  to  dr/dt  exactly  or  a2d(r/rc)/d|  to 
first  order,  the  particle  displacement  can  be  obtained  by  integration.  The 
first-order  particle  displacement  r-r^  or  singly  Ar  associated  with  only  the 
i^h  segment  of  the  signature  is  given  below. 
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(3.62) 


From  this  result  the  particle  displacement  Ar/r0  at  the  beginning  of  only 
the  i^h  segment  (|=2i-2)  is  zero,  of  course,  and  that  at  the  end  of  the  i**1 
segment  (|=2i)  is  given  by  the  following  expression. 


. 2 
tip  r 
^o  o 

3?2I)ir* 


(3.63) 


Hence,  the  total  particle  displacement  Ar/r0  of  the  i^*1  segment,  which  now 
includes  the  displacements  of  all  previous  segments,  is  given  below. 


2 

Ap  r * 
o o 

S^l1* 
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( B^l  Y'1  3/2  r*  a2 
V B+l  ) b+1  ro  aL 
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+ 1 ^ ( 1 - 1 f81-1"0  +l(21-1-s>3)} 


(3*6*4) 


2i-2  < | < 2i 

i = 1,  2,  . . . , oo 

At  large  times  (i  -m»,  £ -»»)  the  motion  of  the  fluid  particle  ceases  and 
the  final  particle  displacement  Ar/r0  is  simply  tsp0T  . For  the 

contact  surface  the  final  displacement  is  Ap0/372Pi*  Note  that  these 
results  ere  in  agreement  with  those  in  Section  3.1  (Eqs.  3.50  and  3*51), 
which  were  derived  by  a different  analysis. 


4.  DISCUSSIONS  AND  CONCLUSIONS 


Acoustic  analyses  have  been  given  for  determining  the  weeik  wave 
motion  in  both  low-pressure-ratio  rectangular  and  pyramidal  shock  tubes. 
Each  sumlysis,  which  was  sufficiently  general  to  deal  with  different  driver 
and  channel  gases  having  different  temperatures,  yielded  a closed-form 
solution  for  the  wave  motion.  Each  solution  was  illustrated  graphically 
with  selected  overpressure  and  particle-velocity  signatures  for  the  wave 
moving  in  the  channel  gas.  Additionally,  solutions  for  the  fluid-particle 
displacement  of  the  channel  gas  and  thus  the  contact-surface  displacement 
were  obtained  for  both  the  rectangular  and  pyramidal  shock  tubes,  and  they 
were  also  shown  graphically  with  selected  results.  Although  the  solutions 
of  the  wave  motion  and  fluid-particle  displacement  for  the  pyramidal  shock 
tube  were  quite  lengthy,  graphical  results  for  the  overpressure,  particle- 
velocity  and  particle  displacement  were  obtained  rather  easily  by  utilizing 
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a relatively  small  digital  computer  (IBM  H30) . Experimental  data  in  the 
form  of  overpressure  measurements,  obtained  in  conjunction  with  suitable 
rectangular  and  pyramidal  shock  tubes,  verified  the  analyses  and  solutions 
for  the  wave  motion  in  both  shock  tubes. 


For  the  special  case  of  identical  driver  and  channel  gases  having 
the  same  temperature,  it  has  been  well  known  that  the  rectangular  and 
pyramidal  shock  tubes  both  produce  a single  finite- duration  wave  or  pulse 
in  the  channel.  For  example,  for  the  rectangular  shock  tube  see  Ref.  l8 
(similar  wave  motion  of  string^,  and  for  the  pyramidal  shock  tube  see  Refs. 

15,  18,  19  and  20.  In  the  case  of  the  rectangular  shock  tube  the  pulse  is 
flat-topped  (see  Fig.  7c),  having  an  amplitude  of  one-half  the  pressure 
difference  across  the  diaphragm  (Ap0/2)  and  a wave  length  of  twice  the 
driver  length  (2rQ) . For  the  pyramidal  shock  tube  the  pulse  is  N shaped 
(see  Fig.  Ilf),  having  the  same  wave  length  of  twice  the  driver  length 
(2r D)  and  an  amplitude  which  is  initially  one-half  the  pressure  difference 
across  the  diaphragm  but  diminishes  inversely  with  increasing  radial  distance 
(Ap0r0/2r) . The  general  solutions  of  the  wave  motion  given  in  Chapter  2 for 
the  rectangular  shock  tube  and  in  Chapter  ^ for  the  pyramidal  shock  tube 
are,  of  course,  in  perfect  agreement  with  these  well-known  results.  However, 
in  the  past  it  has  not  been  fully  recognized  that  in  the  more  general  case 
of  driver  and  channel  gases  having  different  densities,  owing  to  different 
molecular  weights  or  temperatures,  the  amplitude  and  shape  of  the  waves  in 
the  rectangular  and  pyramidal  shock  tubes  can  differ  markedly  from  those  of 
a single  flat- topped  pulse  (see  Fig.  7)  and  N- shaped  wave  (see  Fig.  ll) 
respectively.  In  fact,  this  was  one  of  the  main  reasons  why  the  solutions 
of  the  wave  motion  in  the  rectangular  and  pyramidal  shock  tubes  were  obtained. 

As  explained  in  detail  in  the  Introduction,  a pyramidal  shock  tube 
utilizing  identical  driver  and  channel  gases  having  the  same  temperature  has 
to  be  very  large  (about  200m  long)  if  the  N-wave  produced  in  the  channel  is 
to  simulate  a full-scale  sonic  boom  from  an  SST  aircraft.  Although  a shorter 
pyramidal  shock  tube  could  produce  a similar  duration  simulated  sonic  boom, 
provided  the  sound  speed  of  the  driver  gas  was  relatively  low,  this  simulated 
sonic  boom  would  not  have  the  desired  N shape  and  it  would  also  have  undesirable 
trailing  disturbances.  Consider,  for  example,  the  specific  case  of  sulfur 
hexafluoride  (SFg)  as  the  driver  gas  and  air  in  the  channel.  Owing  to  such 
a low  sound  speed  driver  gas  the  pyramidal  shock  tube  needs  to  be  only  50 % 
as  long  as  one  which  uses  air  in  both  the  driver  and  channel.  However,  the 
distorted  N-wave  which  is  produced  in  the  channel  (see  Fig.  Ilk)  is  quite 
unsatisfactory  for  standard  sonic-boom-research  purposes.  Consequently,  it 
can  be  concluded  that  the  concept  of  using  a low  sound  speed  driver  gas  to 
make  possible  the  construction  of  a relatively  short  and  inexpensive  pyramidal 
shock  tube  for  sonic-boom- research  purposes  is  not  very  feasible. 

At  noted  in  the  Introduction,  a relatively  short  pyramidal-pyramidal 
shock  tube  utilizing  air  in  both  the  driver  and  channel  could  be  used  to 
produce  a simulated  full-scale  sonic  boom  (see  Fig.  5).  This  simulated  sonic 
boom  would  also  not  be  N shaped  and  it  would  also  have  undesirable  trailing 
disturbances  (see  Fig.  3d).  Consequently  the  concept  of  using  such  a short 
and  inexpensive  pyramidal-pyramidal  shock  tube  for  sonic-boom  research 
purposes  is  also  not  very  feasible.  However,  if  a short  and  inexpensive 
sonic-boom  simulator  is  to  be  built,  the  pyramidal-pyramidal  shock  tube  is 
more  practical  than  a pyramidal  shock  tube  which  uses  different  driver  and 


channel  gases.  Firstly,  although  it  can  he  shown  that  the  undesirable 
distortion  in  the  N-wave  produced  by  either  shock  tube  is  about  the  same, 
the  trailing  compression  waves  or  disturbances  for  the  case  of  the  pyramidal- 
pyramidal  shock  tube  (Fig.  3d)  are  much  less  severe  than  the  trailing  steep- 
fronted  disttirbanc.es  for  the  pyramidal  shock  tube  (Fig.  Ilf).  Secondly, 
the  problems  associated  with  the  mixing  of  the  different  driver  and  channel 
gases  in  the  pyramidal  shock  tube  after  the  diaphragm  is  broken  is  nonexistent 
for  the  case  of  the  pyramidal-pyramidal  shock  tube. 

The  analysis  and  solution  for  the  wave  motion  in  a pyramidal- 
pyramidal  shock  tube  for  which  the  driver  and  channel  gases  can  not  only 
be  different  but  also  have  different  temperatures  is  very  similar  to  that 
given  in  Chapter  3 for  the  case  of  the  pyramidal  shock  tube.  For  such  an 
analysis,  if  the  radial  distance  for  the  driver  regior  is  denoted  by  r,  ? 
new  radial  distance,  say  R,  is  needed  for  the  channel.  Therefore  the  length 
of  the  driver  is  rQ  and  the  length  of  the  truncated  portion  of  the  channel 
is  Rq.  Then,  for  example,  the  recurrence  relationship  of  f^(j|)  (see  Eqs. 

3.27  to  3»32)  for  the  pyramidal  shock  tube  can  be  used  for  the  case  of  the 
pyramidal-pyramidal  shock  tube,  provided  the  parameter  l is  modified  slightly 
from  (l+a)/(p£/pi-l)  to  (l+s)/(p2r0/p^R0-l) . Consequently,  the  analysis 
and  solution  for  the  pyramidal-pyramidal  shock  tube  can  be  obtained  quite 
easily  with  this  and  other  slight  modifications  to  the  pyramidal- shock- tube 
results  given  in  Chapter  3.  From  such  a solution  it  can  be  shown  that  the 
overpressure  signature  of  the  wave  in  the  channel  is  generally  not  N shaped 
and  also  has  trailing  disturbances,  very  much  like  those  shown  in  Fig.  11. 

It  is  worthwhile  to  discuss  one  particular  solution  of  the  wave 
motion  in  a pyramidal-pyramidal  shock  tube.  Then  one  can  answer  the 
interesting  question  as  to  whether  a pyramidal-pyramidal  shock  tube  using 
different  driver  and  channel  gases  can  produce  a single  undistorted  N-wave. 
From  the  solution  of  the  wave  motion,  it  can  be  shown  that  the  linear 
variation  in  overpressure  for  an  N-wave  can  occur  only  if  the  parameter  i 
becomes  infinitely  large.  (Note  that  this  condition  is  precisely  that  used 
to  obtain  the  special  solution  given  in  Section  3«1*)  This  condition  on  1 
is  possible  only  if  the  ratio  of  the  driver  and  channel  gas  densities  P2/P1 
is  exactly  equal  to  the  geometrical  ratio  R0/r0.  Consequently,  one  finds 
the  interesting  result  that  the  two  distortions  of  the  N-wave,  one  produced 
by  the  different  driver  and  channel  densities  and  the  other  introduced  by 
the  different  driver  and  channel  divergence  angles,  cancel  exactly  if  P2/P1 
equals  Rg/r0.  Even  though  the  N-wave  is  undistorted,  however,  it  is  followed 
by  a series  of  rapidly  diminishing  amplitude  N-waves  having  the  same  duration 
(2rQ/a2),  exactly  like  the  results  shown  in  Fig.  13.  One  can  conclude, 
therefore,  that  a pyramidal-pyramidal  shook  tube  using  different  driver  and 
channel  gase3  cannot  produce  a single  N-wave. 

There  would  be  no  significant  advantage  for  building  a pyramidal- 
pyramidal  shock  tube  instead  of  a pyramidal  shock  tube  to  produce  an  N-wave 
for  the  simulation  of  a full-scale  sonic  boom.  Firstly,  of  course,  the 
simulated  sonic  boom,  although  it  would  be  N shaped,  woula  have  undesirable 
N-shaped  disturbances . Secondly,  the  overall  length  of  the  pyramidal-pyramidal 
shock  tube  would  not  be  shorter  hut  actually  longer  than  that  for  a correspond- 
ing pyramidal  shock  tube.  It  can  be  shown  that,  for  fixed  diaphragm  and  test 
section  sizes,  the  reduction  in  channel  length  (factor  of  a^/a^)  is  over- 
compensated by  the  increase  in  channel  length,  because  R0/r0  must  equal  Pg/pj^. 
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A technique  first  considered  in  England  (Ref.  9)  was  discussed  in 
the  Introduction.  This  technique  consists  of  using  a low  sound  speed  driver 
gas  in  an  existing  short  pyramidal  shock  tube,  thereby  producing  a simulated 
sonic  boom  with  a longer-than-normal  duration.  Because  the  simulated  sonic 
boom  would  be  an  undesirable  distorted  N-wave  having  trailing  disturbances, 
according  to  the  results  of  Chapter  3,  it  can  be  concluded  that  this  tech- 
nique of  inproving  the  capability  of  an  existing  short  pyramidal  shock  tube 
to  produce  a long  N-wave  is  not  very  practical. 

The  only  type  of  shock  tube  which  produces  a good  N-wave  is,  of 
course,  the  pyramidal  shock  tube,  provided  the  driver  and  channel  gases  are 
identical  and  have  the  same  tenperature.  If  a full-scale  sonic  boom  is  to 
be  simulated  then  one  must  recognize  the  fact  that  a large  pyramidal  shock 
tube  is  required.  It  is  worth  mentioning,  however,  that  smaller  sonic-boom 
simulators  such  as  a loudspeaker-driven  booth  and  a pyramidal  horn  with  a 
mass-flow  valve  can  also  be  used  for  simulation  of  a full-scale  sonic  boom 
(Refs.  3 and  5) • 

The  analyses  and  solutions  of  the  wave  motion  for  both  the  rectan- 
gular and  pyramidal  shock  tubes  is  directly  applicable,  of  course,  to 
finite-sized  planar  and  spherical  Weak  explosions.  This  fundamental  work  is 
relevant  to  the  understanding  of  the  wave  motion  produced  by  such  weak 
explosions,  for  which  the  explosion  gas  is  different  than  the  ambient  gas. 
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TABLE  1 


RECTANGULAR  SHOCK  TUBE 

Franke,  Lursat,  Evrard  Gottlieb  (July  1974) 

and  Devriere  (1972) 


PYRAMIDAL  SHOCK  TUBE 


Webb  and  Pallant  (1968) 

Peter  and  Pfister  (1969) 

Peter  and  Brunner  (1970) 

Slutsky  and  Arnold  (1970) 

Thery,  Peter  and 
Schlosser  (1971) 


Carothers  (1972) 

Ribner,  Morris  and 
Chu  (1973) 

Gottlieb  and  Glass  (1974) 
Gottlieb  (July  1974) 


PYRAMIDAL-RECTANGULAR  SHOCK  TUBE 

Warren  (1966)  Gottlieb  and  Glass  (1974) 

Peter  and  Pfister  (1969)  Gottlieb  (July  1974) 

Peter  and  Brunner  (1970) 


PYRAMIDAL-PYRAMIDAL  SHOCK  TUBE 

Peter  and  Brunner  (1970)  Gottlieb  (July  1974) 

Gottlieb  and  Glass  (1974) 


RECTANGULAR -PYRAMIDAL  SHOCK  TUBE 

Tamboulian  and  Peschke  (1970)  Gottlieb  (July  1974) 

Carothers  (1972) 


SPECIFIC  IMPEDANCE  RATIO  Z (*2P2/a1P1  - Vy2M2/YjMj'  if  Tj  * T2)  FOR  DIFFERENT  DRIVER  AND  CHANNEL  GASES 


00 

cn 

00 

rf 

vO 

CM 

CM 

rH 

sO 

© 

u. 

rM 

rf 

m 

vO 

00 

00 

© 

-a- 

| • 

• 

• 

• 

• 

• 

• 

• 

U 

o 

o 

O 

. 

o 

o 

o 

o 

© 

rH 

10 

to 

i—4 

*-H 

in 

o 

CM 

to 

o 

vO 

U. 

»H 

CM 

in 

nO 

SO 

cn 

o 

rH 

CO 

• 

• 

• 

• 

• 

• 

• 

• 

• 

o 

O 

© 

O 

o 

o 

o 

rH 

rH 

CM 

u. 

rf 

CM 

o 

in 

c- 

o 

Tf 

CM 

rH 

CM 

in 

sO 

c- 

o 

o 

CM 

rH 

• 

• 

• 

• 

• 

• 

• 

• 

u 

o 

O 

o 

o 

o 

o 

rH 

rH 

rH 

u 

CM 

© 

© 

i—4 

i—4 

o 

o 

o 

CM 

O 

H 

CM 

© 

oo 

o 

to 

vO 

CO 

• 

• 

• 

• 

• 

• 

• 

• 

• 

o 

o 

o 

o 

o 

rH 

rH 

rH 

rH 

§ 

CM 

CM 

■*3- 

00 

o 

cn 

m 

vO 

CM 

o 

CM 

to 

00 

o 

rH 

in 

NO 

© 

rJ 

u 

• 

• 

• 

• 

• 

• 

UJ 

O 

o 

o 

tH 

rH 

rH 

rH 

rH 

1 

H 

H 

CM 

00 

to 

o 

to 

00 

CM 

•*> 

c- 

© 

rH 

Tf 

m 

■ 

O 

o 

o 

o 

rH 

rH 

rH 

rH 

h 

n£) 

rH 

o 

00 

© 

rH 

rr 

00 

00 

CM 

Tf 

o 

CM 

rH 

rf 

00 

cn 

CM 

< 

• 

• 

• 

• 

• 

• 

O 

o 

H 

rH 

rH 

rH 

rH 

rH 

CM 

in 

o 

vO 

NO 

to 

00 

to 

C' 

to 

4> 

sO 

o 

© 

tT 

in 

00 

© 

X 

• 

• 

* 

• 

• 

o 

H 

CM 

to 

CM 

to 

rr 

Tf 

in 

o 

O 

C' 

CM 

cn 

rH 

00 

CM 

o 

LO 

00 

00 

in 

to 

cn 

in 

vO 

X 

• 

• 

* 

• 

• 

• 

rH 

rH 

to 

in 

vO 

oo 

u 

2 

O 

to 

o 

to 

00 

© 

to 

00 

© 

o 

© 

c- 

vO 

to 

<d 

• 

• 

• 

• 

• 

*■*> 

cd 

to 

CM 

rH 

© 

o 

o 

O 

o 

©* 

wm 

o 

o 

o 

vO 

to 

vO 

in 

vO 

vO 

to 

CM 

rH 

<n 

in 

• 

• 

• 

• 

• 

• 

• 

© 

© 

fH 

rH 

rH 

rH 

rH 

rH 

rH 

• 

• 

< 

rH 

rH 

o 

oc 

£ 

o 

o 

© 

o 

Tf 

rH 

'tt 

© 

X 

• 

• 

CM 

Tf 

vO 

CM 

© 

HH 

CM 

rH 

rH 

CM 

d. 

O 

1 1 

CM 

U* 

CM 

OO 

H 

M 

CM 

CM 

rH 

10 

Urn  1 

CM 

X 

£ 

3! 

B 

CO 

SO 

u 

u 

SF 

* 

u 

TABLE  3 

PARAMETER  *.  [(1  ♦ z)/(y1z2/y2  -1)  with  Tj  * T2]  FOR  DIFFERENT  DRIVER  AND  CHANNEL  GASES 


200  1.055  0.33  I 0.10  0.14  0.56  0.69  0.82  1.23  3.43  5.54  *od 


FIG 


1 IDEALIZED  OVERPRESSURE  SIGNATURE  OF  A SONIC  BOOM 


c)  PYRAMIDAL-RECTANGULAR  SHOCK  TUBE 


FIG.  2 SHOCK  TUBES  WHICH  HAVE  BEEN  ASSESSED  FOR  SONIC- 
BOOM  SIMULATION  PURPOSES 


FIG.  4 THE  FRANCO-GERMAN  SONIC  BOOM  SIMULATOR 
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ht  1 


FIG.  6 TIME-DISTANCE  DIAGRAM  ILLUSTRATING  THE  WAVE  MOTION 
IN  A LOW-PRESSURE-RATIO  RECTANGULAR  AND  PYRAMIDAL 
SHOCK  TUBE 
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FIG.  7 OVERPRESSURE  SIGNATURES  IN  THE  CHANNEL  OF  A RECTAN- 
GULAR SHOCK  TUBE  FOR -DIFFERENT  DRIVER  AND  CHANNEL 


DISTANCE  (X) 


V. 


FIG. 
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8 MOTION  OF  A FLUID  PARTICLE  IN  THE  CHANNEL  OF  A REC- 
TANGULAR SHOCK  TUBE 


. 1 
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a)  DRIVER  GAS:  AIR 

CHANNEL  GAS:  AIR 

MEASUREMENT  LOCATION:  X > I m 


Ap#  ■ 10  kN/m2 
x#  ■ 46  cm 


4.95 

kN/m2 


b)  DRIVER  GAS:  HELIUM  * Ap  » 10  kN/m2 

CHANNEL  GAS:  AIR 

MEASUREMENT  LOCATION:  X « 2.1  m X.  « 1.5  m 


6.65  /l 
kN/m2 

__ 


k= 10  ms 3-J 


C)  DRIVER  GAS:  AIR  Ap  * 10  kN/m2 

CHANNEL  GAS:  HELIUM 

MEASUREMENT  LOCATION:  X > Im  X.  • 46  cm 


2.8  | 
kN/m2 


* 
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FIG.  9 OVERPRESSURE  MEASUREMENTS  MADE  IN  THE 
CHANNEL  OF  A RECTANGULAR  SHOCK  TUBE. 


DRIVER  GAS:  AIR 

CHANNEL  GAS:  AIR 

MEASUREMENT  LOCATION: 


Ag  » 10  kN/m 


x.  ■ 46  cm 


— PREDICTED  (Z=  1.0) 
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DRIVER  GAS:  HELIUM 

CHANNEL  GAS:  AIR 
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x • 2.1m 


10  kN/m 
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-—MEASURED 
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Ad  ■ 10  kN/r 

CHANNEL  GAS: 
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MEASUREMENT 

LOCATION: 

x * 1 m x»  * 46  cm 
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— PREDICTED  (z  - 2.45) 
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FIG.  10  COMPARISON  OF  PREDICTED  AND  MEASURED  OVERPRESSURE 
SIGNATURES  FOR  A RECTANGULAR  SHOCK  TUBE. 
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FIG.  13  OVERPRESSURE  SIGNATURES  FOR  THE  CHANNEL  OF  A PYRAMIDAL 
SHOCK  TUBE  FOR  WHICH  THE  DIFFERENT  DRIVER  AND  CHANNEL 
GASES  HAVE  EQUIVALENT  DENSITIES  BUT  DIFFERENT  SOUND 
SPEEDS 


PYRAMIDAL  SHOCK  TUBE 


PYRAMIDAL  SHOCK  TUBE 


FIG.  24  PARTICLE-DISPLACEMENT  HISTORIES  AT  DIFFERENT  RADII  IN  THE  CHANNEL  OF  A 
PYRAMIDAL  SHOCK  TUBE 
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GAS : AIR 
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LENGTH  (r.): 


(He) 

Wp.): 

1.58  m 


MEASUREMENT  LOCATION  (r): 


2400  N/m2 
4.6  m 


575 

N/m2 


10  ms 
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DRIVER  GAS:  ARGON  (A) 

CHANNEL  GAS : AIR 

DRIVER  OVERPRESSURE  (Apj,): 
DRIVER  LENGTH  ( r. ) : 1.58  m 

MEASUREMENT  LOCATION  (r): 


2400  N/m2 
7.6  m 
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\ • : 
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TIME  (ms) 


FIG.  27  COMPARISON  OF  MEASURED  AND  PREDICTED  OVERPRESSURE 
SIGNATURES  OF  A PYRAMIDAL  SHOCK  TUBE 


a)  DRIVER  GAS:  CARBON  DIOXIDE  (CO2) 

CHANNEL  GAS : AIR 

DRIVER  OVERPRESSURE  (Afc) : 2400  N/m2 

DRIVER  LENGTH  (r.):  1.58  m 

MEASUREMENT  LOCATION  (r).  7.6  m 


] 


I 
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FIG.  28  COMPARISON  OF  MEASURED  AND  PREDICTED  OVERPRESSURE 
SIGNATURES  OF  A PYRAMIDAL  SHOCK  TUBE 


DRIVER  GAS:  SULFUR  DIOXIDE  (S02) 

CHANNEL  GAS : AIR 

DRIVER  OVERPRESSURE  (Aft):  2400  N/m2 

DRIVER  LENGTH  ( r. ) : 1.58  m 

MEASUREMENT  LOCATION  (r ) : 7.6  m 


FIG.  29  COMPARISON  OF  MEASURED  AND  PREDICTED  OVERPRESSURE 
SIGNATURES  OF  A PYRAMIDAL  SHOCK  TUBE 


DRIVER  GAS:  DICHLOROOIFLUOROMETHANE  (CC12F2) 
CHANNEL  GAS:  AIR 

DRIVER  OVERPRESSURE  (Aft,):  2400  N/m2 

DRIVER  LENGTH  (r.):  1.56  m 

MEASUREMENT  LOCATION  (r):  7.6  m 


i 


FIG.  30  COMPARISON  OF  MEASURED  AND  PREDICTED  OVERPRESSURE 
SIGNATURES  OF  A PYRAMIDAL  SHOCK  TUBE 


a)  DRIVER  GAS:  OCTAFLUOROCYCLOBUTANE  (C4F8) 

CHANNEL  GAS:  AIR 

DRIVER  OVERPRESSURE  (Ap#):  2400  N/m2 

DRIVER  LENGTH  (r.):  1.58  m 

MEASUREMENT  LOCATION  (r)  : 7.6  m 


FIG.  31  COMPARISON  OF  MEASURED  AND  PREDICTED  OVERPRESSURE 
SIGNATURES  OF  A PYRAMIDAL  SHOCK  TUBE 


APPENDIX  A 


,th 


The  following  expression  for  the  i velocity  potential  f.'(Tj)  was 
derived  in  Chapter  3 (Eq.  3.29).  1 


A i 

fi(t|)  = a2P2(a+l)  I 

3=1 


i-l 

3-1 


g-1 

B+l 


i-3 


2 

H+l 


3-1 


I’(ri-2i+l) 


If  Ap0r 2l/ £i2P2(H+l) j (b-1)/(b+1)  and  2/(h+1)  are  denoted  by  A,  B and  C 
respectively,  the  first  five  velocity  potentials  can  be  expressed  as  follows. 


f{(r,)  = A[  I^(t)-l)] 

f'(ri)  = A[B  I£(t)-3)  + C I2(t)-3)] 

f^(ri)  = A[B2I^(ti-5)  + 2BC  I2(tj-5)  + C2I^(tj-5)] 

f£(ri)  = A[B3I£(t)-7)  + 3B2C  I2(tj-7)  + 3HC2I'3(t\-7)  + C2I^-7)] 

f^(ri)  = A[BUq(n-9)  + 4b3C  I2(t|-9)  + 6b2C2I3(t)-9)  + 4bc3I,;(ti-9)  + C4I’(tj-9)] 

The  3 function  Il(tj-2i+l)  derived  in  Chapter  3 (Eq.  3-32)  can  be  expressed  in 
the  following  form. 


Jp(w)  = [(l+i)(u)/£)  /2 : + a(u>/£)]exp(-ui/i) 


J^(w)  = [(l+£)(uj/i)3/3:  + a(<V4)2/2'.]exp(-u>/i) 
J£(w)  = [(l^)(Vi)V;  + «(u/i)3/3’.  ]exp(-aj/i) 


APPENDIX  B 


The  following  expression  for  the  overpressure  in  the  channel  of  a 
pyramidal  shock  tube  utilizing  different  driver  and  channel  gases  was  derived 
in  Chapter  3 (Eqs.  3.37,  3*38  and  3«39)« 

00 

i=l 


gjd) 


S(sxan*y 

J=1 


Jj(i-2i4€) 


i = 0 


i = 2,  3, 


• • • 9 


00 


l + (l+i)exp(-u?/i) 


1+i  f U) 


— \J)  - TJ=TTT  V I 
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J-l- 


J = 0 


exp(-u>/J)  j = 1,  2, 


If  (b-1)/(b+1)  and  2/(b+1)  are  denoted  by  B and  C respectively,  the  first  five 
functions  for  g” ( t ) can  be  expressed  as  follows. 


g£(0 

" Jo<*> 

g£(0 

= C J”d-2) 

ggU) 

= BC 

+ c2j^(i-4) 

Slid) 

= b2c  j£(s-6) 

+ 2BC2j^(|-6) 

+ c3j^U-6) 

g^U) 

= B3C  j£( 1-8) 

+ 3B2C2Jg(|-8) 

+ 3bc3j^(i-8)  + c4j£U-8) 

The  first  five  functions  of  Jj'(?-2i+2)  are  listed  below. 

Jq(w)  = - l + (l+i)exp(-u/i) 

J^(«)  a [(!+/)(<*)//)  + B]exp(-w//) 
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